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ABSTRACT 


The effect of the adhesive layer on crack propagation 
In composite materials is Investigated. The composite me- 
dium consists of parallel load carrying laminates and buffer 
strips arranged periodically and bonded with thin adhesive 
layers. The system Is assumed to approximate boron-epoxy 
composites. The strips, which are assumed to be Isotropic 
and linearly elastic, contain symmetric cracks of arbitrary 
lengths located normal to the Interfaces. Two problems are 
solved for both plane strain and plane stress cases. In 
the first problem, thin adhesive layers are approximated by 
uncoupled tension and shear springs distributed along the 
Interfaces of the strips for which only the case of Internal 
cracks can be treated rigorously. The second problem Is 
Introduced In order to study the case of broken laminates 
and to detect the true singular behavior In the presence of 
the adhesive layer. In this case the adhesive Is treated 
as an Isotropic, linearly elastic continuum. General ex- 
pressions for field quantities are obtained In terms of In- 
finite Fourier Integrals. These expressions, with relevant 
boundary and continuity conditions, give a system of sing- 
ular Integral equations 1n terms of the crack surface 
displacement derivatives. By using appropriate quadrature 
formulas, the Integral equations reduce to a system of 
linear algebraic equations which Is solved numerically. 


Thent stress Intensity factors and some significant stress 
components are calculated. The results are compared to 
those obtained by neglecting the adhesive layer (perfect 
bonding assumption). Results are presented In graphical 
form. 
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I. INTRODUCTION 

Fracture mechanics is the study of the strength of a 
structural member that contains a crack. A normally ductile 
member may behave in a orittle manner if it contains cracks 
or other flaws which are suff iciefi tly large. Adhesive- 
bonded composites tend to have flaws due to the complexity 
of shape* chemical dissimilarities, and assembly procedures. 
These flaws, under load, may develop into cracks with a re- 
sulting brittle failure. In applying fracture mechanics, 
it is assumed that all real structures have initial flaws 
or cracks, and that failure is caused by the propagation of 
the largest of these. The techniques of fracture mechanics 
can measure the intrinsic toughness of the material, which 
determines the load-carrying capacity of the structure in 
the presence of flaws. Therefore, if the size of the 
largest flaw in a particular structure is known, minimum 
toughness standards can be established for the materials in 
the structure. One has to adopt a proper fracture criterion 
and decide on the type of "load factor" to be evaluated. 
Usually the fracture criterion consists of a simple com- 
parison between a calculated load factor and a material 
constant which is determined from certain standard experi- 
ments. 

Inmost fracture analyses, two basic (essentially equiv- 
alent) approaches have been used with variations. With one 
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approach, due originally to Griffith (1920), the energy 
required to propagate a crack of a given size Is considered. 
A crack will propagate If the rate of release of the stored 
energy per unit growth of the crack exceeds the rate of 
change of the surface energy required by the new surfaces. 
According to Irwin [1], the stress field In the vicinity of 
a crack tip can be adequately defined for studies of crack 
extension by a single parameter, proportional to the stress 
Intensity factor. Since this parameter Is a function of the 
applied load and crack size (it increases with load), when 
the Intensity of the local tensile stresses at the crack 
tip attains a critical value, a previously stationary or 
slow-moving crack propagates rapidly. This critical value 
defines the "fracture toughness" and It Is a constant for 
a particular material, since cracking always occurs at a 
given value of local stress intensity regardless of the 
structure In which the material has been used. Fracture 
toughness has the same relationship to brittle design that 
yield strength has to ductile design. Variations of these 
theories also have been useful. When a significant degree 
of plasticity takes place In the structural member, the 
usefulness of the elastic stress intensity factor as a 
correlating parameter becomes questionable. For cases where 
large-scale yielding can be expected, there are other cor- 
relating parameters that have been suggested In recent 
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years such as the J*1ntegral» crack opening displacement, 
plastic stress and strain Intensity factors, etc. 

In the fatigue of bulky structures with no high stress 
concentrations, the major portion of the fatigue life Is 
elapsed before the formation of macrocracks. The remain- 
ing portion of the fatigue life Is relatively very short. 

On the other hand, 1n composite structures the formation of 
a dominant macrocrack may take placi relatively early In 
the fatigue life. Hence In such cases the propagation 
phase, I.e., the number of load cycles necessary for the 
fatigue crack to reach a critical length at which the struc- 
ture may fall statically, represents the major portion of 
the total fatigue life. Therefore In composite structures 
propagation and arrest of fatigue cracks Is a major subject. 
Since the stress Intensity factor Is the simplest and the 
most appropriate single variable used 1n studying the fa- 
tigue crack propagation. Its evaluation attracts consider- 
able attention. 

Because of the ever-increasing use of modern composite 
materials In a wide variety of structural applications, the 
mechanics of multi-phased materials has attracted consider- 
able attention, particularly within the last ten or fifteen 
years. A great amount of work in this area has been con- 
cerned with the influence of the localized imperfections 
on the overall response of the medium regarding its failure. 
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The primary Interest is mostly In the initiation and propa- 
gation of fracture in the composite material. The fracture 
process may start as the initiation of a fatigue crack at 
a local imperfection in the composite material. This crack 
then propagates with the cyclic effect of applied loads re- 
sulting in structural failure at stress levels considered 
moderate in relation to the theoretical strength of the ma- 
terial. There are two types of failure: (1) due to a con- 

trolled rate of cracking, arising from a steady rate cf 
stressing; (2) due to catastrophically fast crack growth. 

The latter one needs more attention since the growth cannot 
be controlled easily. The main reason for the use of rel- 
atively low stiffness and high toughness buffer strips par- 
allel to main load-carrying laminates in designing with 
high strM',^!th composite materials, is to improve the fatigue 
crack ;/ropagation and arrest characteristics of the struc- 
ture. 


The use of very strong epoxy type adhesives has been 
very common in joints of flight and space vehicle structures 
in which lightweight and high fatigue strength are dominant 
requirements. Epoxy based adhesives are also being used 
increasingly in stiffening, joining, and repairing precast 
prestressed concrete and other structures. 

The adhesive layers, which serve as the bonding agent 
in composites, have not been treated adequately in 
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literature (see [2]). In the past some of the problems» 
particularly the problems relating to the traction-free 
boundaries, have presented considerable analytical diffi- 
culties, There are some finite element solutions which are 
good to the extent that they are reliable. However, these 
solutions miss the correct singular behavior near the cor- 
ners and hence, from the viewpoint of fracture studies, 
they are of limited value. In recent years, with the in- 
troduction of the concept of generalized Cauchy kernel and 
the development of the related numerical techniques (see 
[3], [4]) it seems that some of these problems can now be 
treated in an analytically correct manner. For example, 
the problem of a composite isotropic plate which consists 
of oarallel load-carrying laminates and buffer strips has 
been solved by Erdogan and Bakioglu [5]. In this study 
the effect of the adhesive layers has not been taken into 
account, Eisenmann and Kaminski [6] had considered this 
problem before. They concluded that crack arrest could be 
achieved through the use of buffer strips in the primary 
load-carrying laminates. However, their analytical work 
in evaluating the quantities which are useful in design 
considerations is not complete. The same problem has re- 
cently been considered by Delale and Erdogan [7] for ortho- 
tropic materials again neglecting the adhesive layers. 
Erdogan and Civelek [8] have treated the thin adhesive 
layer as a shear spring in the contact problem for a thin 
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elastic reinforcement bonded to an elastic plate. In the 
stress analysis of a metal ba *e plate stiffened by a fiber- 
reinforced composite layer Erdogan and Arin [9] assumed 
that the two materials are bonded through an adhesive which 
Is treated as a two-dimensional shear spring. The problem 
of an elastic plate bonded to a rigid horizontal substrate 
through an adhesive layer has been considered by Williams 
[ 103 . In this work, the adhesive layer has been assumed to 
react as a Winkler elastic foundation. It should be noted 
that the essential assumption introduced by Winkler allows 
for vertical motion and dilatation stress only. In one of 
the most recent works, Updike [11] has Investigated the 
effect of adhesive layer elasticity on debonding of a blis- 
ter test specimen. The adhesive layer has been treated as 
a distributed spring or elastic foundation which transmits 
normal and shear stresses between plate and support. 

The objective of this work is to Investigate the effect 
of adhesive layer on crack propagation In composite ma- 
terials. A composite medium, which Is generated by parallel 
main load-carrying laminates and buffer strips bonded 
through thin adhesive layers, Is considered. Materials of 
main laminates and buffer strips are assumed to be Iso- 
tropic and linearly elastic. Main laminates and buffer 
strips are arranged periodically and they are assumed to 
contain symmetric fatigue cracks of arbitrary lengths and 
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traction-free surfaces normal to their longitudinal direc- 
tion. The composite medium Is loaded In y-directlon con- 
siderably far from crack region (see Figure 1). The solu- 
tion to this problem can be obtained by superposition of 
solutions for the following two problems; (1) A strip 
having no crack loaded In y-dlrectlon. and (2) a strip 
having a crack whose surface Is subjected to the negative 
of the stress distribution obtained at the same location In 
the first problem resulting from the applied loads. Solu- 
tion of the first problem Is relatively simple and straight 
forward hence one pays more attention to the second problem 
(see Figure 2). Therefore we solve the perturbation prob- 
lem in which crack surfaces are subjected to prescribed 
tractions. 

Thin adhesive layers are approximated by distributed 
uncoupled tension and shear springs. As it can be seen 
In the relevant references mentioned above, short of con- 
sidering the adhesive as an elastic continuum, this Is the 
most sophisticated model for adhesives in literature. A 
formulation is given for both plane strain and plane stress 
cases. General expressions for displacement and stress 
components are obtained by solving field equations using 
Fourier transform technique. Applying boundary conditions 
and the continuity conditions at the interfaces a system 
of singular Integral equations in terms of crack surface 
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displacement derivatives Is derived. By using appropriate 
quadrature formulas, these Integral equations are converted 
to a system of linear algebraic equations which Is solved 
numerically. Stress Intensity factors and stress components 
for Imbedded cracks are computed and results are given In 
Figures 3-16. The results are compared to those obtained 
without taking thickness of the adhesive layer Into account. 
Then the case of broken laminates Is considered. It Is ob- 
served that the spring model approximation Is not suitable 
under these circumstances when the crack touches Interface. 
Therefore. In order to be able to examine the singular be- 
havior of the cleavage stress, we Introduced the problem 
described In Chapter III. In this problem, the adhesive Is 
treated as an Isotropic and linearly elastic medium. All 
laminates are of the sam'^ material and thickness but the 
cracks are of different lengths (see Figure 17). The prob- 
lem Is solved for Imbedded cracks In order to determine 
limitations for the spring model approximation In the first 
problem. Comparison of two solutions can be seen In Fig- 
ures 18-23. The case of broken laminates can now satisfac- 
torily be solved. The power of singularity at the crack 
tip is determined from the characteristic equation obtained 
by following Muskhel Ishvll 1 [12]. The Integral equation Is 
replaced again by a system of linear algebraic equations 
and this system Is solved to calculate the stress Intensity 
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factor and stresses. These results are given in Figures 
24-37. Discussion of numerical results and conclusions 
can be found in Sections 2.5, 3.6. and Chapter IV* 
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II. SPRING MODEL PROBLEM 


Z.'t Formulation of the Problem 

Consider a medium consisting of Infinitely many number 
of linearly elastic and Isotropic strips bonded together. 
These strips constitute periodically arranged load carrying 
laminates and buffer strips which are of different thick- 
nesses and material properties. There are fatigue cracks 
In laminates and buffer strips normal to the Interfaces. 

Main laminates and buffer strips are bonded through thin 
layers of adhesive. Loads are applied away from the crack 
region and In a direction parallel to the strips {see Fig- 
ure 1). Solution for a strip having a crack of traction- 
free surface and loaded sufficiently far from the crack 
region can be obtained by superposing the solutions of 
(1) a loaded and uncracked strip, and (11) a strip having 
a crack whose surface Is subjected to the negative of the 
stress distribution acting on the Image of the crack in 
the uncracked strip due to remotely applied loads. There- 
fore In solving the abovementi oned elastostatic plane prob- 
lem, one should first solve the problem of a strip having 
a crack whose surface is subjected to a prescribed traction. 
In this work we will consider the singular part of the 
problem In which the self-eqjll Ibrating crack surface trac- 
tions are the only external loads (see Figure 2). 
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First, we will establish the forms of the field quan- 
tities for a strip having a crack. Forms of the exp esslons 
for the field quantities for a strip having a crack and 
which Is loaded sufficiently far from the crack region can 
be obtained by superposing the general expressions of these 
quantities for: (1) an uncracked strip, and; (2) an Infin- 

itely large medium having a crack. 

The basic equations for a linearly elastic. Isotropic 
medium In plane problems can be written as (see [13]): 


Strain-displacement relations: 


e 


XX 


du 

ax ’ 


e = 
yy ay ’ 



M 

ay ax * 


(2.1a-c) 


where e , Yyw are component." of the strain tensor, 

AA yy Ay 

u, V are displacements in x- and y-directlons In a Cartesian 
coordinate system. 

Stress-displacement relations: 


1 _ _ (c+l au . 3 -k av 

- 2'(K-r) w 3y ■ 



3-IC au . K+1 iv 

zXk-T) ax 2 (k-i ) ay 


j_ 

2p’^xy 


1 ;ljy. + 

2 'ay 3 x' ’ 


(2.2a-c) 
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where a^y, are components of the stress tensor, p 
Is the shear modulus, ic>3-4v for plane strain case* and 
tc*(3-v)/0+v) for generalized plane stress case, v being 
the Poisson's ratio. 

Equilibrium equations, in the absence of body forces: 

u3ii 4. 2. 3 _ n 

(2.3a,b) 

2 ^ /3U ^ ivw ^2 Q 

ay 'ax ay' ’ 

where ^ ^ . 

2.1.1 Havana no Cn,ack 

Taking the Fourier cosine transform of Eq. (2.3a) and 
sine transform of Eq. (2.3b) in y-direction, and combining 
the resulting equations, we obtain: 



2s2 


d^v 

BP- 


+ s"9 


= 0 . 


(2.4) 


where (') stands for the sine transform and s is the trans- 
form variable. Note that the strip is symmetric about both 
X- and y-axes. The solution of Eq. (2,4) is: 

v= [f(s) +--^ g(s)]s“^cosh(sx) + g(s)x sinh(sx) , (2,5a) 


where f(s) and g(s) are unknown functions. Similarly, 

u = - [f(s) g(s)]s"^sinh(sx)-g(s)x cosh(sx) , (2.5b) 
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where ('^) Implies the cosine transform. Taking the Inverse 
transform of Eqs. (2.5), displacement components are found 
to be 

u(x,y) ■ - |•|^^[f(s)- g(s)]s"^s1nh(sx) 

+ g(s)x cosh(sx)}cos(sy)ds , 

- (2.6a,b) 

v(x.y) * |- Jq([^(s)+ ^ g(s)]s"^cosh{sx) 

+ g(s)x s1nh{sx)}s1n(sy)ds . 


Substituting Eqs. (2.6) Into Eqs. (2.2) one obtains the 
following expressions for the stress components: 


= - Ij [f(s)cosh(sx)+ g(s)sx s1nh(sx)]cos(sy)ds , 


^<Jyy(x,y) |^([f(s)+ 2g(s)]cosh(sx) 

+ g(s)sx sinh(sx)}cos(sy)ds , 


(2.7a-c) 


«00 

^Tjjy(x.y) s I ^{[f(s)+ g(s)]s1nh(sx)+g(s)sx cosh(sx)}s1n(sy)ds . 


2.1.2 Ela6t^c. Solid Having a. CKack 

Taking now the Fourier sine transform of Eq. (2.3a) 
and cosine transform of Eq. (2.3b) In x-direction and com 
b^n'!ng again the resulting equations, one obtains 
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d'^a d*0 . . n 

37^ - 2r‘ gyr + r**Q » 0 , 


( 2 . 8 ) 


In which r is the transform variable. Note again that the 
medium Is symmetric about both x> and y.*axes. Solution to 
Eq. (2.8) vanishing when becomes then 

a * - m(r)r’^(^ - ry)e"*^^ , (2.9a) 

where m(r) Is an unknown function. Similarly, 

V - m(r)r”^i^ + ry)e"*'^ . (2.9b) 

Taking the Inverse transforms of Eqs, (2.9), we get 

u(x,y) * - ~ [ m(r)r’^(-*^ - ry)e"*^^sin(rx)dr , 

° (2,10a,b) 

v(x.y) - I rm(r)r'’(Sjl + ry)c'*'^co5(r*)dr . 

Substituting Eqs. (2.10) Into (2.2) one obtains the stress 
components as follows: 

2ll ^xx^^*y^ * ■ TT J^«i(r){l-ry)e’’"^cos{rx)dr , 

^ cfyy(><»y) * - “ |^m(r)0+ry)e“’*^cos(rx)dr ,{2.11a-c) 

^ Tjjy(x.y) « - I j^m{r)rye’’"^s1n{rx)dr . 

2.1.3 Su}(3&n.poiltion 

Now consider a strip which has a crack on which there 
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are prescribed tractions as the only external loads. Gen- 
eral forms of dlspUcement and stress components are es- 
tablished by superposing those for an uncracked strip* and 
a cracked space as mentioned before. Hence we obtain: 

u^(x^.y) ■ - I |^{s'^[f^(s) g^(s)]s1nh(x^s) 


x^g|(s)cosh(x^s)}cos(ys)ds 
- f |*r“^m^(r)(-^ - ry)e"'''ysin(x^r)dr , 
v^(x^.y) « I j^{s"^[f^(s) + g^(s)]cosh(x^s) 

+ x^g^(s)s1nh(x^s))s1n{ys)ds 

tfiff ^ 

+ I * <f*y)e'*^^cos{x^r)dr , 


(2.12a, b) 


as the general expressions for displacement components for 
the i-th strip where l^l for the main load carrying lamin- 
ates and 1=2 for buffer strips. Similarly, general ex- 
pressions for the stresses become: 


y^*^1xx^^r^^ ' - |^[f^(s)cosh(x^s)+sx^g^(s)s1nh{x^s)]cos(ys)ds 

- fOO 

-|J m^(r)(l-ry)e’’''^cos(x^r)dr , 
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( 


• I |^{[f^(s)+2g^(s)]cosh(x^s) 

* sx^g^(s)s1nh(x^s)}cos(ys)ds 

- ff j^n»^(r)(l+ry)e‘*^^cof5(x^r)dr , 

(2.13a-c) 

57 " 0[fl(s)*9i(s)]5inl>(x,s) 

•f sx^g^{s)cosh(x^s)}s1n(ys)ds 

- ^ |^(n^(r)rye"'*^s1n(x^r)dr* . 


The unknown functions f^(s), g^(s) and m^{r), (1*1,2) can 
be determined by using appropriate boundary and continuity 
conditions. 

2.1.4 C ontinuity ConcUt^onA 

Consider the elastostatic plane problem shown In Fig- 
ure 2. Load carrying laminates (^) having thickness of 
2h^ are bonded to buffer strips of thickness 2hg by 

means of thin adhesive layers (^) of thickness 2h2. Main 
laminates and buffer strips contain symmetrical cracks of 
lengths 2a and 2b respectively located on the line y=0, 
normal to y-axis. Thin adhesive layers are approximated by 
distributed tension and shear springs. That is, the x 
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component of normal stresses and the x-y component of shear 
stresses In the adhesive layer are assumed to be constant 
through the thickness 2h3. Then the continuity conditions 
may be written as: 


0, ,,(!>, .y) ■ 02 x*(-" 2>1') • 0<y<- 

T,xy(h,.y) ■ Tj^y(-hj.y) , 0 <y<- 


(a.l 4 a,b) 


Ujdijty) - UjC-hj.y) • - r 

,1° ( 2 .) 5 a,b) 

v,(h,,y) - Vjt-hj.y) ■ - ^ ■ 0<y<“ 


where hQ«2h3, Wq“^ 3* generalized plane 

stress case Ce3yy“0, 

(I-2V3) for plane strain case (£3yy*0t ^322*0) where E Is 
the Young's modulus and the subscript (3) stands for the 
adhesive layer. 

Substituting Eqs. ( 2 . 12 ), { 2 . 13 ) into Eqs. ( 2 , 14 ), 
( 2 . 15 ), taking Fourier sine and cosine transforms of the 
resulting equations In y, and making use of the Integral 
formulas given by (A.l) in Appendix A we obtain the follow- 
ing four equations: 

K,-l 

( +2\^shQa2 ) f ] ( s )+«3f2 ( s )+( ^‘l +sh-j 012 + 2 ^ h^h^ ) g^ ( s ) 

K«— 1 

+ ot3+sh20i^)92(^) “ ” Fi 1 ( s) “2 X^shQF| 3(5 ) -F 21 (s ) , 
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K,+l 

(a2+2X3Shjja^)f^(s)-a^f2(s)+(-^ a2+sh|a^+2X3S*hQh^02+2X3Shpa^)g^(s) 

K *^1 

- (-^ a4+sh2a3)92(s) ■ -F^2^s)‘*'2^3®^o ^14^® ^*‘’^22^®^ * 

-a2f^ (s)+A2a^f2^®^’®^'*l9l (s)+^2®^392^®^ “ Fi3(s)-X2l^23^s) • 

(2.16s-d) 

{s)+X20t3f2(s)+(ot^+sh^ot2)9^ (s) 

+ X2(a3+sh2a^)92(s) ■ Fi4(s)+X2F24(s) , 

where 

» sinh(sh^) , 
ttg » cosh(shi) , 

* sinhtshg) t 
= cosh{sh2) » 

and 

X^ “ y^/P2 * X 2 “ * 

X 3 = . X^ = Pi/Eq . 

The functions Fjj(s), (1=1 t2; j=1-4) are defined by (B.l) In 
Appendix B. Solving Eqs. (2.16) simultaneously we obtain 
the unknown functions f^(s), g^(s), (1=1,2) In terms of 
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Infinite Integrals of the other unknown functions m^(r), 
(1=1,2) as follows: 

= I [a..F, ^(s)+b. ,F2.(s)]e"®^^ , 

(2.17a.b) 

9i(s) » I Cc^jF^ j(s)+d^jF2J(s)]e■®^^ , (1-1.2) . 

J “ 1 

Known functions a^j, b^j, c^j, d^j ( 1 = 1 ,2;j = l -4) are given 
In Appendix B by (B. 2) . 


2.1,5 BoundaK^ Cond^t.ion6 

The boundary conditions related to symmetry about x- 
and y-axes have already been used in deriving the general 
expressions for displacement and stress components. The 
remaining boundary conditions are used to determine the two 
unknown functions m^ ( r), ( 1=1 ,2) . These boundary conditions 
may be stated as follows: 

(2.18a,b) 

V2{X2»0) = 0 , b<|x2l<h2 

<^ 1 yy ( X-j ,0) = ’•P'|(Xl) , 0 <|x^j<a 

(2.19a.b) 

O2yy(x2.0) = -P2(X2) . 0<|x2l<b 


Eqs. (2.18) can be replaced by: 


^ v^(x^,0) = 0 . a^<jx.|<h^ . (a^=a,a2=b) 


( 2 . 20 ) 
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and 




p -§|- v^(x^,0)dxi - 0 . (1»1.2) . (2.21) 

“"*1 

The unknown functions ')» (1s1,2) are determined 
from the mixed boundary conditions (2.19) and (2.20). 

2 . 2 Derivation of Integral Equations 

Substituting Eqs. (2.12b), (2.13b), and (2.17) Into 
Eqs. (2.19) and (2.20) we could obtain two sets of dual In- 
tegral equations for m^(r), (1*1,2). However, we will fol- 
low a more direct procedure to solve the problem by defining 
new unknown functions, G^(x^), as 

^l^’^l^ “ v^(x^.O) , (1 = 1,2) ; 0<jx^|<h^ . (2.22) 

From Eqs. (2.12b) and (2.22), 

Kj+1 , fOO 

Gl(Xi) = (- f J^m^(r)s1n(x^r)dr . (2.23) 

Taking the Fourier sine transform of Eq. (2.23) In x^ and 
using Eq. (2.20), 

mi(r) = - (t)s1n(rt)dt (2.24) 

is obtained. 

Now using Eq. (2,24) and the integral formulas given 
by (A. 2) in Appendix A, the functions defined by (B.1) may 
be expressed as: 
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> 




(2.25) 


where 




and 


l-K. 


Ml 


(s,t) = -s(h^-t) + , 


1+K. 


N^2(s.t) = s(h^-t) - 
N^ 3 (s,t) “ s(h^-t) , 

N^ 4 (s.t) = s(h^-t) - 1 , (i=l,2) 


(2.26) 


(2.27a-d) 


From Eq. (2.13b) we can write, 


^ ~ |^{[f^(s)+2g^(s)3cosh(x^s)+sx^g^(s)s1nh(x^s)}ds 

- Hm [“ f m^(r)n+ry)e"*"^cos(x^r)dr] . ( 2 . 28 ) 

y + 0 + ^ •'0 ' ^ 

With m^(r) as defined in Eq. (2.24), the second integral in 
Eq. (2.28) becomes; 


( 1 = 1 . 2 ) . 


(2.29) 


N^te that, 


- 23 - 


due to symmetry of v^{x^,y) about y-axls. 


Substituting Eqs. (2.17), (2.25), and (2.29) Into Eq. 
(2.28) and using Eq. (2.30) 

J"' ' -a 

+ c^jSX^slnh(x^s)] 

u 4)i. rb /“ 

J^M2.(s.t)[(b^j+2d^.)cosh(x^s) 


+ d^jsx^s1nh(x^s)]e‘ + -^dt,(1=1,2) 


(2.31) 


Is obtained. Now Eqs. (2.19) with Eq. (2.31) give, 


If"®!'*! .. . f" 


? J t-)T, 

-a * 


dt + 


1+<I rt> 


-a 


^11 (^1 >t)Gi { t)dt + jif “ J k*j 2(x*| »t)G2(t)dt = 

2 -b 


1+K, 




p,(x,) . -a<x,<a , 


1 

n 


j-b 02^^^ J ^■*■'^2 


-b 


t-x.- 


dt + 


1+K 


1 '-a 

1+K 


(2.32a,b) 

fb 

l<2i(x2,t)G^(t)dt + J I<22(x2tt)62(t)dt = 


^2^^2^ * "b<X2^b 


as the Integral equations for the new unknown functions 
( t) , ( i =1 , 2 ) where the kernels k^j(x^,t), (1,J=1,2) are 
given by 




(2.33) 


ds . 

The expressions for the functions K|j(x^,t,s) are given In 
Appendix B by (B.5). 

The dominant parts of the Integral equations (2.32) 
have a simple Cauchy type singularity for a<h^ and b<h 2 . 
However, It Is worthwhile to take a closer look at the 
Fredholm type kernels k^j(x^.t). The Integrands of the 
kernels vanish as s-*-« and are bounded everywhere except for 
s=0. From (B.5) It can be seen that as s->-0 the functions 
L^j(x.j,t,s) = K^j(x^ ,t,s)e"®^’^‘^"^^ behave as s"^. That Is, 
for smal 1 val ues of s , 

L^j(x^,t,s) = ^ + e^(l) , (2.34) 

where i,j=l,2, and are constants Independent of t. We 
can write, 

.00 -e .00 

L^j(x^ ,t,s)ds * J j(x^ ,t,s)ds + J j (x^ ,t,s)ds 

(2.35) 

for any 0<e«”. Choosing a very small e, and hence replacing 
L^j(x^.,t,s) In [0,e] by Its asymptotic expansion around s = 0 
[see Eq. (2.34)], we can rewrite Eq. (2.35) as 

j(x^ ,t,s)ds = 0'(1 )]ds + I j(x^ ,t,s)ds . 

(2.36) 
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1 


I 


Now consider the last equation and» for example* the fol- 
lowing Integral appearing In Eq. (2.32a): 

f k, , (X|,t)G, ( t)dt . |“g, (t)dt{ 9(1 ) ]<ls 

+ I K^^(x^,t.s)e"®^^^"^^ds} 

£ 

= I E|jjO(l)ds + |^K^^(x^.t,s)e’®^'^^"^^ds]G^(t)dt 
-a 

+ j%^(t)dt • ^ds. (2.37) 

-a 

Using the single-val uedness condition. Eq. (2. 21). we can 
see that the last Integral In Eq. (2.37) vanishes. The re- 
maining Integrals are all bounded and they can be evaluated 
numerically. One can show similarly that the singularity 
around s=0 Is removable for all kernels. 

All the known special cases may be recovered from 
Eqs, (2.32) by letting a.b.h^hg. and h^ (or E^.p^) go to 
proper limits. For example, for Eqs. (2.32) become 

Identical to those found in [5]. Note that making h^ very 
small has the same effect as making E^, \iq very large due 
to spring model assumption. That is, a very small h^ cor- 
responds to a very stiff adhesive and as h^^ increases the 
r.dhesive becomes less stiff. Hence the case of direct ad- 
hesion of laminates to buffer strips [5] can be achieved 


- 26 - 


by letting Infinitely large will then 

Imply that there Is almost no bond between laminates and 
buffer strips. Therefore, letting hg-^* (or Eq.Uq-^O) we re- 
cover the case of a center-notched strip having traction- 
free straight boundaries [14]. By letfing h^-^O (or Eq.Uq-^“’) 
and making laminates and buffer strips be Identical (a»b, 
hi'hg. obtain the solution for the case 

of collnear cracks In a homogeneous medium [15]. If we let 
b*0, h^'O (or and hg-*-®, the Integral equations 

(2.32) reduce to an Integral equation In 6^(t) for the prob- 
lem of two half planes bonded through a center-notched 
strip [16]. 

2 . 3 Solution of Integral Equations 


2.3.1 Mo CA.a.c.k In Baj^zfL Stfi^pA [a<h^ ,b^0) 

In this case the integral equations (2.32) 
, fa G,(t) fa 1+K, 

? J dt + J k^ ^ (x^ ,t)6i (t)dt = - (x^ ) , 

-a ‘ -a ' 


reduce to: 

-a<x^<a 

(2,38) 


where Gi(t) Is unknown, and Pi(x-|) and k^i(xi,t) are known 
functions which are Holder-continuous In the closed interval 
[-a, a]. A function f(x) is said to be Holder-continuous in 
[-a, a], if for any two points x-j .Xg in [-a, a] the following 
condition is satisfied; 
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Iftxg) - f{x^)| < Ajxg-x^l^^ , -a<x<a 


(2.39) 


I 


where A and y are positive constants, and 0<y<1. 


6^(t) and k^^(x^,t) being Hd1der>cont1nuous , the sec- 
ond Integral In Eq. (2.38) gives a bounded function of x^ . 
Hence, the singular behavior of G^(t) may be obtained by 
studying the dominant part of Eq. (2.38) only: 

1 fa Gij(t) 

?J (2.40) 


where F^(x^) contains all the bounded terms In Eq. (2.38). 

Define the sectlonally holomorphic function, 

1 fa Q,(t) 

Assume that the unknown function G^(t) has the following 
form: 


6^(t) = 4>^(t)(t+a)^(-tfa)“ 


• s 

(t+a) near t = -a , 
(-t+a)“ near t = a , 


(2.42) 


where 4>-|(t) Is a bounded continuous function in the closed 
interval [-a, a], and a, 6 are yet unknown constants restrict- 
ed by -l<a,e<0, which implies that the unknown function 
G.|(t) has an integrable singularity at the end points t =ta. 

Following [12] and using Plemelj formulas [17] given 

below 
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(2.43) 


(X]) “ 'I'l (j<1 ) ® » 

X -1 6,(t) 

'^l(><,) + '^^(x,) •TTl, 


the value of (x^ ) around » .a and x^ ■ a can be found as 

(j>^(-a)(x^+a)^+4>^(-a)i;/^^(x^)+i{'^2^’‘i) around x]=-a, 


1 cot 7T( 


(2.44) 


(2a) <|,^{a)^^a)%^(a)i|)^ 3 (x^)+i|)^ 4 (x^) around x^-a 


where 'i'^^(x^), (i»l-4) are bounded In [-a. a]. F»*om Eqs, 
(2.40), (2.41). and (2.44) 


-(2a)“cot7r3 (|i^ (-a)(x-|+a)®+(2a)^cot ira(()^ (a)(-x^+a)“= , -a<x^<a , 

(2.45) 


where Pi](^l) contains all the bounded terms. Multiplying 
Eq. (2.45) first by {-x^+a)“® and letting x^-+a, and then by 
(x.j+a)’^ and letting x^-^-a 


cot na » 0 
cot tt3 » 0 


(2.46a,b) 


are obtained for -l<a,3<0 since ( -a ) , c|)^ ( a ) ?* 0. 

Eqs. (2.46) give 

a = 0 = -l/2 . (2.47) 


The index of the problem which is defined by 


I 


K » - (a+3) , 


(2.48) 


Is related to the physical nature of the problem. In this 
problem, has Integrable singularities at both ends 

and K«+l. From Eqs. (2.42) and (2.47) we obtain 

■^l ( t ) 

6i(t) » . -a<t<a . (2.49) 

I /a'-t=" 

The solution of the problem will contain one arbitrary con- 
stant [12]. Theoretically this constant Is determined by 
using the single-val uedness condition, Eq. (2.21). However, 
due to symmetry, Eq. (2.21) will be automatically satisfied 
if one considers a solution satisfying Eq. (2.30). 


Now substituting Eq. (2.49) Into Eq. (2.38) we obtain 


1^3 4>i (t) 1 l+K-j 

J-J -(2.50) 


Define non-dimensional variables aj,T by 


X, - au , t = ax , -a<x,,t<a 


(2.51 ) 


Then Eq. (2.50) takes the form 

I I ^5= Ci + irak,^(au),aT)]dT 
IT J.| /PF ^ 


- 1 , -1 <W<1 , 


(2.52) 


where 


4y ij).(aT) 

= T^- • -br 


(2.53) 
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dnd P^(xi) Is assumed to be (x<| )« « constant. Actually* 

this Is not an assumption. The sel f >equ1 1 Ibrating pressure 
acting on the crack surface Is p^ * constant when the strips 
are subjected to uniformly distributed normal loads In y- 
dlrectlon. The Integral In Eq. (2.52) can be evaluated by 
using the Gauss-Chebyshev Integration formula given In [17], 
Hence, Eq. (2.52) becomes 
N , 

I 6|(t^)[ "^ - :— ' + TTak^^(aaij,aT^)] = -N , (2.54) 

1 ® 1 1 J 

where 

» cos[(2i-l )tt/2N] , (1 = 1 N) . 

(2.55a,b) 

Wj » cos(Jtt/N) , (j=»l , . . . ,N-1 ) , 

are the roots of related Chebyshev polynomials, Eq, (2.54) 
provides (N-1) linear algebraic equations for (N) unknowns, 
0 ](t^), (1 = 1,. ..,N), The single-valuedness condition, Eq. 
(2.21), which can now be written as 
N 

1 6 (t,) = 0 (2.56) 

1 = 1 ' ’ 

completes the system of (N) equations for (N) unknowns. 
However, one can take advantage of symmetry and can deal 
with less number of equations. Using Eq. (2.30) we can 
write Eq. (2.52) in the following form: 
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I 




"1 • 


(2.57) 


where 


ic^^U.x) » a[kn(auj,aT) - kuCaai.-at)] . 


Ml 


11 


(2.58) 


Now Eq. (2.57) can be replaced by 

0 1.1 


• -N • <J*’ "> • 


(2.59) 


In which and wj are the same as In Eqs, (2.55) and 2n»N. 
Eq. (2.59) constitutes a system of (n) linear algebraic 
equations in (n) unknowns, Q^(x^), (i = l,...,n). The single^ 
valuedness condition, Eq. (2.21), Is automatically satis- 
fied by the solution obtained from Eq. (2.59). 

After determining 0-|(x) at discrete colocation points, 
the field quantities can be computed numerically. In frac- 
ture problems, one Is Interested mostly In the so-called 
"stress Intensity factor". The stress intensity factor 
may be defined 1n terms of cleavage stress, Oiyy(x<|,0), 
and may be expressed in terms of 0^(x) i\s follows (see 
Appendix C) : 

= lim a^yy(x^.O) = - p^0^(l) , a<h^ . (2.60) 
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I 


2.3.2 Ho CAack -in Hain Laminat&i {a<’0,b<k2) 

The Integral equations (2.32) reduce to 

1 6M(t) 1+K« 

7 j -feu <“ ♦ J l<22(«2.t)G2(t)dt. - P2(«2) . -Il<*2<b . 

4. y y 


-b 


(2.61) 


Following a procedure similar to the one followed In sec- 
tion 2.3.1, we determine the behavior of the unknown func- 
tion G 2 (t) as 


G„(t) 
2 


-b<t<b 


(2.62) 


where <j»o(t) Is a bounded continuous function in the closed 
interval [-b,b]. 

Substituting Eq. (2.62) into Eq. (2.61) one obtains 


, fb 4>o(t) 1 1+K, 

’’ • ^2.63) 


Defining non-dimensional variables ntt by 

Xg * bn . t = bT , -b<X 2 ,t<b (2.64) 

Eq. (2.63) takes the following form 

If [^ + ^bk««(bn,bT)1dT = -1 , -l<n<l . (2.65) 

where 
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( 2 . 66 ) 




4Mg 4-2(bT) 
T+^ . ^ 


and p 2 (x 2 ) ■* P 2 ' constant. Using the Gauss-Chebyshev Into 
gratlon formula given In [17], Eq. (2.66) becomes 





+ 7rbk22(bnj .bt^ )] 


m 


-N . 


(2.67) 


Note that 


'j = “j 


(j»l N-1) 


( 2 . 68 ) 


T^,o)j being defined In Eqs. (2.55), and hence Eq. (2.67) 
can be written as 




1 


■ 0) , 


+ TTbk22(UWj ibT^ )] = -N 


(2.69) 


The singl e-val uedness condition, Eq. (2.21), In this case, 
becomes 

1 = 1 ^ 


Eqs. (2.69) and (2.70) constitute again a system of (N) 
linear algebraic equations for (N) unknowns, 02^^1^’ 

(1=1,. ..,N). Taking again the advantage of symmetry, 
namely using Eq. (2.30), Eqs. (2.69) and (2.70) can be re- 
placed by 

J ’''‘22*“^ ’ ’ 

1=1 1 J 1 J (2.71) 


\ 
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where 


= b[k22(^>“j»bT^) - k22(bWj,-bT^)3 , (2,72) 

end 2n » N. Eq. (2.71) give (n) linear algebraic equations 
for (n) unknowns, C>2^'^i^* (1“l.**‘.n), ^2 obtained from 

Eq. (2.71) will satisfy the single-val uedness condition, 

Eq, (2.21), automatically. 

The stress Intensity factor, in this case, may be ex- 
pressed as (see again Appendix C): 

kjj = 11m /2Tx^^ a 2 yy(x 2 , 0 ) = - P 2 e 2 (l) , b<h 2 . (2.73) 

X2*^b 


2.3.3 CA.ack'i Lam^natai and SuiUe/i. ia<lij,b<h 2 ) 


The singular behaviors of G^(t) and G 2 (t) at the end 
points t = -a,a and t = -b,b respectively will be the same 
as in the cases where there are cracks in one type of strips 
only. Now substituting Eqs. (2.49) and (2.62) into Eqs, 
(2.32) we obtain 


1 fa ‘<»i(t) 1 1+K, rb (})2(t) 


1+Krt /-a (|), (t) , fu ^2'“' 1 

ia k2^(x2,t)dt + ~ -^= 


" Pl(X]) . -a<x^<a , 

, ,b ♦g**’ , 1 


(2.74a.b) 


1+K, 




Using the non-dimensional variables defined by Eqs. (2.51) 
and (2.64). Eqs. (2.74) and (2.21) become 

" f! 7T^ (ei(T)[Jjj*mkll(a«i.aT)] + 92(T) ^ X,nbk,j(ai..bT)}dT 


-1 , -l<tlJ<l t 


(2.75a.b) 


» f.Ti^ X2«kj,(bn.aT) t ej(T)[:^+itbk22(bn.bT)])<lT 


« - 1, -l<n<l t 


J ej^(T)dT = 0 , (k=l ,2) . 


(2.76) 


Use of Gauss-Chebyshev Integration formula for the evalua- 
tion of Integrals In Eqs. (2.75) and (2.76) l.eads us to 
N , Po 

1^1 ^0 * 

(2.77a,b) 

N P 

^®1 V^4l < •3'^i J ^ ^ ^bl<22(ba)j ,bT^ ) ]> = -N , 

(a=i N-1) 


I 0.(t.) = 0 . (k-1,2) 
i=l ’ 


(2,78) 


which constitute a system of (2N) linear algebraic equations 




for (2N) unknowns, 0 |^(t^), (k=l,2;1 = 1 N). We can use 

Eq. (2.30) with Eqs. (2.75) to obtain 

(2.79a,b) 
(j=l .... .n) 


where 

kigCujjT^) ® b[k^ 2^®***j >1 ) ” k^ 2(®^j )3 * 
i^2‘] (^j ) ~ a[k2i (bWj ,ax^ ) - k2i(bWj,“aT^)] . 


(2.80a,b) 


Eqs. (2.79) give an (NxN) system. 

Definitions of the stress intensity factors will be 
the same as those in Eqs. (2.60) and (2.73) since the s1ng< 
ular behaviors of the functions 6^(t), (1=1,2), at the end 
points are not affected by the cracks in the other strips 
for a<h^ and b<h 2 - 

In these numerical solutions we cannot determine the 
values of the unknown functions 0]»02 points di- 

rectly from the system of linear equations, since we are 
not allowed to choose the end points as colocation points, 
as it can easily be seen in Eq. (2.55a). Nevertheless, 
using an appropriate summation formula [18], which is based 
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on the properties of Chebyshev polynomials, we can deter- 
mine 0 ^( 1 ), 62 ( 1 ) In terms of the values of 6 i *02 
location points. From [18], we can write 

1 n s1n[^J^ (21-1 )ir] 

V^1> • <'‘=’•2) • (2-81) 

One should keep In mind that the above approximate formula 
can be used for a= 0 =-l /2 only. 


2.4 Stresses at the Interfaces 


We consider the stresses at the Interfaces for the 
case In which there are cracks In main laminates (a<h-j,b= 0 ) 
only. The other cases can be treated similarly. Since the 
shear stresses, and the normal stress In x-directlon, 

o^xx* assumed to be the same for both strips at the 
Interfaces, for the sake of simplicity we consider the 
stresses In the buffer strips when there are cracks In the 
main laminates. 


Now, from Eqs. (2.13) we can write 


02xj^(-h2.y) = - I j [f2(s)cosh(h2s) + sh2g2(s)s1nh(h2s)]cos(ys)ds , 


2pr 


®2yy^"*^2’^^ "I * 2g2(s)]cosh(h2s) + sh2g2(s)s1nh(h2s)} 

cos(ys)ds , 
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1 


i T2xyl-"2-y) 


/QO 

I J^{[f2^s) + g2(s)]sinh(h2s) + sh292(s)cosh(h2S)} 

sin(ys)ds , {2.82a-c) 


since m 2 (r) is zero by definition* Eq. (2.24), for b=0. 
Using Eqs. (2.17) with Eq. (2.25) 


f 2 (s) * I e"®^2|^[M^j(s,t) -M^j(s,-t)]6^(t)dt , 

j •"! 1 

^2^^^ ° ^ e"^^2|^[M^j(s,t) -M^j(s,-t)]G^(t)dt , 

j “1 1 


(2.83a.b) 


are obtained. Then substituting Eqs. (2.83) into Eqs. 
(2.82), the stress expressions become 


2I2 ^2xx^"^2*y^ “ ~ in%T * 

?J^‘^2yy^-^2»y^=¥n%r ’ 

^T2xy(-h2.y) = - |^[h23(y.t) - h23(y,-t)]G^ (t)dt , 

(2.84a-c) 


where 

h2|^(y.t) = |^H2|^(y,t.s)e"®^^T^^ds , (k=l-3) 


(2.85) 


and 


“Shr 


i2i(y.t,s)= 5; Nij(s,t)[a2jC0sh(h2s) + C2jSh2Sinh(h2s)]e cos(ys) 
J ™ i 
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H22^y»^*s) ® I N^j(s,t)[(a2j+2c2j)cosh(h2s)+C2jSh2Sinh(h2s)]e"®^2 

cos(ys) , 

H23(yit,s) =^^^N^j(s,t)[(a2j+C2j)s1nh(h2s)+C2jSh2Cosh(h2s)]e’^^2 

sin(ys) . (2.86a-c) 

Now using Eqs. (2.49) and (2.53) with the non-dlmen- 
slonal variables ud.t defined by Eq. (2.51), and defining p 
by 


y = h^p , (2.87) 


we obtain 



1 f 

°2xx^"*’2»^^/Pl *“‘^2 ? J( 

1 6i(t) . , , 


1 0 

c^2yy^-f’2»y)/Pr^2?J, 

6i(t) 

fi22(P«'^)d'r . 

(2.88a-c) 

T2xy(-h2»y)/Pi = -^2 ¥ 1 

' h23(p.x)dT . 



where 


= a[h2k(biPiaT) - h2k(h^P,-aT)] , (k=l-3) . (2.89) 


The Gauss-Chebyshev Integration formula given in [17] again 
can be used to evaluate the integrals in Eqs. (2.88). Hence 

Xg n 

‘^2xx^"^2’^^^Pl " ” IT 01 ^'^1^*^21 ’ 
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(2.90a-c) 


02yy(-h2.y)/Pi = T ’ 

T2^(-h2.y)/p, - - f Jj9,(T,)ii23(p,T,) . 

2.5 Numerical Results 

The analysis carried out is valid for both plane 
strain and generalized plane stress cases. The systems of 
equations (2.59), (2.71), and (2.79) are solved numeric- 
ally for 9 |^(t^.), (k»l,2;i = l n). In computing the 

kernels, the infinite integrals are evaluated by using the 
approximate Laguerntquadrature formula [19] since the in- 
tegrands have exponential ly decaying behavior. Once ®|^(”'^j) 
are determined, the stress intensity factors and the 
stresses at the interfaces can then be computed easily. 

Note that the stresses are computed for the perturbation 
problem whose configuration and loading conditions are 
shown in Figure 2. Assuming that there is no constraint 
in x-direction, and the composite medium shown in Figure 1 
is loaded in y-direction sufficiently far from the crack 
region (i.e., the dimension of the medium in y-direction is 
large compared to that in x-direction), the crack surface 
tractions in Figure 2 satisfy the following conditions: 

p,/p- = E,/E« for plane stress case, 

’ ^ ' (2.91a,b) 

P 1 /P 2 = ( 1 -v^)/E 2 ( 1 -v^) for plane strain case. 
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I 


Here ly and E 2 are the Young's moduli for the strips ^ 
and respectively. The problem is solved for three 
combinations of materials: 

Combination I: 

Vi-j * 6. 65^2 • * '^l” ^*3 ~ ® ’ ^2 ~ 0.45 . 

= 12 . 38x10^ Vm2 (17.955x 10^0 . 

E 2 = 2 . 0 xl 0 ^°N/m^ (2.9xlO®psD . 

£ 3 = 0.31 xlO^Vm^ {0.45xl0®psi) . 

Combination II: 


^2 ” V^i » V*2 * ® • 025u^ I ~ ^2 ” ^2 ~ ^ ' 

E^ « £2 = 12.38xl0^%m^ (17.955 x 10®psi) , 

£ 3 = 0.31 xlO^°N/m^ (0.45xl0®psi) . 

Combination III: 

^1 “ ^2 f ~ ® ^ ^^^2 * ^1 "" ^2 • V3“0.35 . 

E^ =£2 = 2.0 xlO^°N/m^ (2.9xl0®psi) . 

£3= 0.31 xlO^°N/m^ (0.45xloSsi) . 

Among these. Combination I is the most significant 
one which is assumed to approximate boron-epoxy composites 


4 
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f 


having buffer strips of the same material with different 
stiffness. In this combination, the first strip Is the 
stiffest one and the adhesive layer Is the softest one. 
Figures 3-16 show the results obtained for material Combin- 
ation I. Figure 3 shows the variation of the normalized 
stress Intensity factor, p^ , with a/h^ for 

hg'h^' b = 0 (section 2.3.1), h 3 /h-| = .00, .05, and .10; 
for plane strain and generalized plane stress cases. It 
may be observed that: (a) Increases with Increasing 

h 3 /hi ; (b) K.j Is larger In generalized plane stress case; 
(c) Increases with Increasing a/h^. Figure A Is the 
plot of Kg = kjjZ/B' pg vs. b/hp for h^ = hg, a = 0 (section 
2.3.2), hj/hg = .00, .05, and .10; for plane strain and ^:en- 
erallzed plane stress cases again. One can conclude that; 
(a) Kg increases as hg/hg Increases; (b) Kg Is larger again 
In generalized plane stress case; (c) The trend In varia- 
tion of Kg with b/hg depends on the ratio h 3 /hg. For small 
values of h 3 /hg (up to '-.04) Kg decreases with increasing 
b/hg. For larger hg/hg ratios Kg increases with Increasing 
b/hg. It is well known that the stress intensity factor at 
the tip of a crack approaching the Interface of two differ- 
ent materials increases if the crack is in the stiffer 

material and decreases if the crack is in the softer 
Tin 

Problem is formulated in such a way that the ratio h 2 /hi 
can vary. But the results are shown for h 2 /h] = 1 to con- 
form with the experimental programs carried out elsewhere. 
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material. When hg/hg is small, it seems as if the crack in 
buffer strip approaches the main laminate, which is stiffen 
than buffer strip, when b/hg increases. But when hj/hg is 
large, the adhesive layer, which is softer than buffer 
strip, seems to dominate. One should also keep in mind 


that increase in h^/hg is equivalent to relative decrease 


in £ 3 , M 3 due to spring model approximation. 


Figures 5-8 


show the variation of and Kg with a/h^ and b/hg for 


hg s hi , h 3 /h^ = .00, .05, and .10. Main laminates and buffer 
strips both contain cracks (section 2.3.3). We can observe 


and state similar conclusions as in Figures 3 and 4. In 


Figures 9 and 10, variation of K-j with h 3 /h^ is shown for 
hg = h-j, a/h-|“0.8, b/h^= 0 . 0 , 0 . 8 , and h 3 /h^ = 0-100. As 
h 3 /h^ increases, also increases tending to the asymp- 
totic value for h 3 '>“> which corresponds to the case of 
center-notched infinite strip having traction-free straight 
boundaries, since h 3 -*-“> is equivalent to E 3 ,M 3 -^ 0 . Figures 
11 and 12 show similarly the variation of Kg with hg/hg 
for h^ = hg, a/hg = 0 . 0 , 0 . 8 , b/hg = 0 . 8 , and h 3 /hg = 0 - 100 . 

K^ and Kg become equal to each other when h 3 -+« (or E 3 ,M 3 *^ 0 ). 
Note that the solution (stress intensity factor) for a 
center-notched single strip is independent of the material. 
All the results obtained for h 3=0 are exactly the same as 
those in [5], 


r 


- 44 - 


In Figures 13,14 distributions of the stress compon 


ents f^2yy ^2xy Xg » -hg boundary for generalized 

plane stress case with hg = h-| , h 3 = 0.05h^, b = 0 are shown. 
As the crack propagates, the stresses at the Interface In- 
crease. Note that «»ixx^^1 " ®2xx^’^2’^^* ^Ixy^^l’^^ ’ 

Tgxy(-hg,y) and <^xx’%y symmetric whereas Is anti- 
symmetric about the x-axIs. Figures 15 and 16 show the 

variations of ^2yy^"*’2'®^ '^2xy^"*'2*^1 ^ ^S^in 

for generalized plane stress case in which hg = h^ , b»0. 
Beginning at h 3 = 0, Og^^ first Increases (up to hg^.OUh^) 
and then decreases as h 3 /h<| Increases, limiting value being 
zero for h 3 ^“, whereas Xg^^^ decreases continuously, vanish- 
ing for h 3 -*-«. This Is expected since h 3 -><» (or equivalently 
relieves the constraints on the boundaries. 


We attempted to solve the problem for a * h-j , b = 0. 

But some unavoidable troubles arose. For example, a singu- 
larity power of -1/2 Is found from a characteristic equa- 
tion which does not contain any material constants. It Is 
now obvious that the spring model used to approximate the 
adhesive layers is not suitable for this case. The model 
is valid for cases where crack tip 1s away from the ad- 
hesive layer. So we introduced the problem described in 
the next chapter in order to be able to examine the effect 
of the adhesive on the solution of the problem when a “ h-j . 
In this problem, the adhesive is treated as an elastic 
continuum. 
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III. CONTINUUM MODEL PROBLEM 


3.1 Formulation of the Problem 

Consider the following problem: Infinite number of 

strips of the same material and the same thickness of 2h^ 
are bonded through adhesive layers of thickness The 

main laminates contain periodical ly arranged symmetric 
cracks of lengths 2a and 2b on the line y =0 perpendicular 
to y-directlon (see Figure 17). The strips ^ and 
are symmetric about x- and y-axes. Therefore the general 
expressions derived In Section 2.1.3 for displacement and 
stress components are still valid for the main laminates 
In this problem. But the adhesive layers are symmetric 
with respect to x-axIs only. Therefore we have to derive 
displacement and stress expressions for the adhesive layer. 

3.1.1 Viiplac.&ment and StfLz.i4> EKpfie.&6ion6 th& 

AdhzAivz Layzfi 

Now let us consider Eq. (2.4). The adhesive layer Is 
symmetric about x-axIs only. The solution to Eq. (2.4) for 
a strip with no symmetry about y-axis becomes 

V3 = ^ f4(s) + 94(5) + 8X393(5)] s1nh(sx3) 

(3.1a) 

+ ~ [f3(s) +-^ 93(5) + sx3f4(s)] cosh(sx3) , 
where f ^ ( s ) ,g^ ( s ) , ( 1 = 3 ,4) are unknown functions. We can get 
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U 3 • - [^3(5) - 


( 3 , 1 b) 


K--1 

93(5) + sx3f^(s)] sinh (SX3) 


1 ^3 * 

- j [94(5) --g— f4(s) + 5X393(5)] 00511(5x3) * 


in a similar way. Taking the inverse transforms of Eqs. 
( 3 . 1 ), the displacement components are found as, 

U3(x3,y) • - 1 1 {[f3(s) - 93(5) + sx3f^(s)] sinh(sx3) 

+ [94(5) -—g— f4(s) + 5x393(5)] cosh(sx3)}cos(sy)ds , 

( 3 . 2 a,b) 

V3(x3*y) = 1 ? ([“V* 5X393(5)] sinh(sx3) 

K *^1 

+ [f3(s) +-|— 93(5) + sx3f4(s)] cosh(sx3)}sin(sy)ds . 
Substituting Eqs. ( 3 . 2 ) into Eqs. ( 2 . 2 ) we obtain, 

fCO 

= - Iy[f 3 ( 5 ) + sx3f4(s)] cosh(sX3) 

+ [94(5) + 5x393(5)] sinh(sx3)}cos(sy)ds , 


57: ®3yy<*3'^> = I | ([fsts) + 293(5) + sXjf4(s)] cosh(sX3) 

+ [94(5) +2f^(s) + 5x393(5)] 5inh(5X3)}cos(sy)ds , 
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r3j^y(x3.y) - 1 * 93 ^^^ * sx3f^(s)]sinh(sx3) 

( 3. 3a-c) 

+ [f 4 (s) + 94 ( 5 ) + sx393(s)]cosh{sX3)}s1n(sy)ds 

for a strip having no crack. 

3.1.2 BoundaA.^ and Cont^nuitij Conditions 

Boundary conditions on the line y = 0 for this problem 
are the same as those In Eqs. (2.18). (2.19). However, the 
continuity conditions do change since the adhesive layer is 
not modeled as distributed springs any more. Taking the 
adhesive layer as an elastic continuum we can write the 
following conditions of continuity: 

oixx^^ivy) = «?3xx^'^3*y^ * 

Tixy(hi.y) = X3xy(-h3.y) . 

u^(h^.y) = U3(-h3,y) , 
v^(h^,y) = V3(-h3,y) , 

(3.4a-h) 

®2xx^"h*y^ = ’ 

T2xy(‘h^»y) “ * 

V 2 (-h^,y) = V 3 (h 3 *y) , 0 <y<«> . 
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Substituting Eqs. (2.12)» (2.13)» (1”1*2)» and Eqs. 
(3.2)» (3.3) Into Eqs. (3.4); then taking the Fourier sine 
and cosine transforms of the resulting equations In y; and 
using the Integral formulas given by (A.l) In Appendix A we 
obtain the following eight equations: 

\a2f^(s)+Xsh^ot^g^(s)-a^f3(s)-sh3a3g3(s)+sh3a^f4(s)+a3g^(s) « -XF^3 , 

Xa^ f^(s)+X(a^+sh^a2^9](s)+a3f3(s)+(a3+sh3a^)g3(s)-(a^+sh3a3)f^(s) 
-<* 494 ( 5 ) * XF ^4 , 

<*Tf^(s)+C(l-ic^)<*l/2+sh^a2]g^(s)+a3f3(s)+[(l-K3)a3/2+sh3a4]g3(s) 

-Csh3a3+(l-K3)a4/2]f4(s)-a^g4(s) « -F„ . 

<*2^1^s)+[sh^a^+(l+o^)a2/2]g^(3)-a^f3(s)-[sh3a3+(l+K3)a4/2]g3(s) 

+[(l+K3)ot3/2+sh3a4]f4(s)+a3g^{s) = -F^2 • 

Xa2f2(s)+Xsh^a^g2(s)-a4f3(s)-sh3a3g3(s)-sh3a^f4(s)-a3g^(s) = -XFgg , 
Xa^f2(s)+X(a^+shia2)g2(s)+a3f3(s)+(a3+sh3a4)g3(s)+(sh3a3+a4)f4(s) 

aif2(s)+[0-Ki)£*^/2+sh^ag]g2(s)+a3f3(s)+[(l-K:3)a3/2+sh3a^lg3(s) 

+[sh3a3+(l-K3)a4/2]f4(s)+a4g^(s) = ♦ 

<*2f2(s)+[sh^a^+(1+ic^)a2/2]g2(s)-a^f3(s)-[sh3a3+(l+K3)a^/2]g3(s) 
-[(l+K3)<*3/2+sh3a4]f4(s)-a3g4(s) = -Fgg . 
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(3.5a-h) 



where 


a-] = slnh(sh^) , 
ttg cosh(sh^) , 
aj * s1nh(sh3) » 

- coshish^) , 

and X = u^/P 3 , The functions F^j(s), ( i = 1 .2 ; j = l -4) are 
defined in Appendix D. After solving Eqs. (3.5) simultan- 
eously, the unknwon functions f^(s), g^(s), (i=l-4) are ex- 
pressed in terms of infinite integrals of m^(r), {i=l,2) 
as follows: 

-s hk 

f<(s) = I ta,jF,j(s) + b,jF2j(s)]e 

(3.6a,b) 

g^(s) j(s) + d^jF2j(s)]e ^ , (i = l-4) . 

where k = l for i = l,2 and k *■ ) for 1 = 3,4. a^j, b^j, c^j, 
d^j, (i,j=l-4) are given in Appendix D. The unknown func- 
tions m^(r), (1=1,2) can again be determined from the 
mixed boundary conditions (2.19) and (2.20). 

3 . 2 Derivation of Integral Equations 

Defining G^(x^), (1=1,2) as in Eq. (2.22) we can simi- 
larly obtain 
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I 


m^{r) = - UGi(t)sin(rt)dt , (1=1,2) . (3.7) 

Following a precedure similar to that followed In section 
2.2 and using the pertinent equations one can write 

|^M,j(s.t)[(a,j«c,j)cosh(x,s) 

-Shl 1 /b r» 

+ c^jsx^s1nh(x^s)]e Gg(t)dt J M^j(s,t) 


'Shl, 


[(b^j+2d^j)cosh(x^s) + d^^sx^s1nh(x^s)3e ds) 
4y^ 1 r^i Qjtt) 

K,+l v\ t-x. • 

1 ia,^ 1 


(3.8) 


Now substituting Eq. (3.8) Into stress boundary conditions, 
Eqs. (2.19) . we obtain 


i" J J " "^jir Pl^^l 

-a I ^•b I 


) . 


-a<x^<a 


I k2i(x2,t)G-j(t)dt + | j [^>^2 ^ h ^’^l ^ 

-a "h 

(3.9a,b) 


1+K, 


•b<X2<b 


where the kernels k,jj(x,.,t), (i,j = l,2) are given by 


k^j(x^,t) = I K^j(x^,t,s)e 


-s(hi-t) 


ds , 


(3.10) 
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I 


and K^j(x^,t,s) are given In Appendix D by (D.5). k^j(x^,t) 

behave similar to those used In Section 2 . 2 . 

Again some special cases can be recovered from Eqs. 
(3.9) by letting some length and/or material parameters go 
to proper limits. For example* for a>b one can obtain the 
one-crack case in [5]. For a = b and hg = 0 or , 

V 3 = we recover the problem of colinear cracks in an elas- 
tic solid [15]. If we let p^-^O the integral equations (3.9) 
reduce to an integral equation for the problem of a center- 
notched strip [14]. By making h^*^” we obtain the case of 
two half planes bonded by a center-notched strip [16]. 
Finally by letting h.j->«> or we recover the case of a 

crack in an infinite elastic solid. 


3 . 3 Solution of Integral Equations 


Following exactly the same procedure as in Section 
2.3 we can obtain, for the case of crack in strip ^ 
only (a<h^ ,b= 0 ) , 

. (j=l.....n) (2.59) 


which is to be solved for 0 -|(t^), (i = l,...,n) and 9 i(t) is 
defined in Eq. (2.53). For the case of crack in strip 
only (a» 0 ,b<h.|), we will have again 
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I 


I 92<^i)CrV* 

1=1 ‘ “j 




+ nk 


22 


(«j.T^)3 


-M . (j«l....,n) (2.71) 


in which 


02^^) 


4y^ <}>2(bT) 

1+K.j * 


(3.11) 


In case of crack in both strip ^ and strip @ (a<h^,b<h^), 

"'‘ll(“j-^l)] * 92 (t,)it|;, 2 (( 0 j.t,)) = -N. 

(j=l.....n) (3.12a,b) 


are the equations to be solved for B|^(x^), (k=l ,2;i = l , . . . ,n)i 
Eqs. (2.60) and (2.73) are still valid for the stress In- 
tensity factors for a<h-| , b<h-| In this problem. 

3,4 Case of Broken Laminates (a=h-j »b=0) 

3,4.1 The. Xnte.QKaZ Equation 

In this case, Eqs. (3.9) reduce to 

tF J +TTkil(x^,t)]G^(t)dt = - p^(x-,) , -h^<x^<h^ . 

-h] 

(3.13) 

However, k^-|(x^,t) is no longer bounded in -hi<Xi ,t<h^ and 
it has point singularities at t = h-j and Xi = +h^ . Eq. (3,13) 
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Is not an ordinary singular Integral equation since It 
contains generalized Cauchy kernel. Integrand of the in- 
finite Integral giving k^^(xpt) Is bounded and continuous 
everywhere In 0<s<«, and the singularity near s = 0 is re- 
movable (see Section 2.2). The divergence in the Integral 
is due to the behavior of the integrand as One can 

separate the asymptotic part of K.|i(x^,t,s) by writing 

^ (xi I t»s ) = ] s ^ K-|i^(Xtt»s) I (3.14) 

where 


K^ 1 s(^l»tiS) “ 11m K^‘|(X‘^tttS) I 


(3.15) 


and K^^f(x^,t,s) is bounded in 0<s<“. 

k^^g(x^.t) = |^K^^g(x^,t,s)e ^^ds 

r“ ^ -s(h,-t) 

k^^^(Xift) “ J 1 f ^ ^ 


Now define 


(3.16a,b) 


so that 


kii(x^^it) ~ kii^(x.|it) 


(3.17) 


k^l^(xpt) being bounded in -h^<Xi,t<h|. After taking 
limits it Is found that 

7rK^ls(xpt,s) = {[4Q^h^(h^-t)s2-Q.|(8h^-6t)s + 3Q^ +Q2lcosh(x^s) 

- 2Q^x^[2(h^-t)s^-s]sinh(Xis)}e , (3.18) 
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where 


Q, = (1-X)/(X+Ki) . 

Q 2 * “Xk2)/{ 1+Akj) • 

Using the following integral formula [20] 


(3.19) 



X|S. 


} ds 


d” r 


(3.20) 


we can evaluate the integral in Eq. (3.16a) and obtain 
irk^ls(xi.t)=]- [4Q^(hi+Xi)" ^ + 12 Q^(h^+x^) 3 ^+ 30^ - Qg] 

^t^(2h-j+‘x^ 

+ 7 [4Q,(ll,-x,)" 3 ^*^'’i'‘52^ 


ttq?).',— ,1^ • 

Now substituting Eq. (3.17) into Eq. (3.13) 

u f +^k,, (xpt)]G^(t)dt+ \^^^(x^ ,t)G^(t)dt 

-hi 1 ihi 


(3.21 ) 


1+ic 


- p^(x^) , -h^<x^<h^ 


(3.22) 


is obtained. The terms in the bracket constitute a typical 
generalized Cauchy kernel. 
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1 


3.4.2 Ckan.aatzfLlAtic Bciuatio n 


Assume that G-|(t) has an Integrable singularity at 
ts+h^ which can be expressed as 

^ttIy 

lt|<h^ (3.23) 




((>^(t)e 


*1 


(hj-t )Y (t-h,)Y(t+h,)Y 


where 0<Re(Y)<l and <J)^(t) Is Hdlder-contlnuous In the 
closed Interval -h^<t<h^. In order to determine the power 
of singularity, y» one should study Eq. (3.22) near t = +h^ , 
Consider the following sectlonally holomorphic function 
rl^l Gi(t) 


, ;in u, VI 
-hi 


dt 


(3.24) 


which becomes 


(h'j ) 


4*1 ( ”hi ) 1 ■ I • I ■ I /I 

^ ° (ah,)^ ^ (2+h,)'' ' (Zh^ji^sImiY (z-h,)'*' 


(3,25) 


by using Eq. (3.23) and following [17]. Here 4*^o(z) is 
bounded everywhere except at the end points +h^ where 

A 


14*10(2)1 


i^thiTl ’ 


Re(Yi)<Re(y) 


(3.26) 


A is a real constant. From Eq. (3.25), using the Plemelj 
formulas given by Eq. (2.43) we obtain 


- 56 - 


1 


1 


where the behavior of Is similar to that of We 

can similarly evaluate 


~ f ^G^(t)(h^tx^) — ~ — I— 

TT ih^ 1 11 dx^ t Uh^+x^J (2h^)Ys1mtY (h^+x^)^ 





dt = . JiLvilL .. 

(2h^)Ys1n7TY (h^+x^)^ 


+ 






(3.28a,b) 


where the behavior of i/i^g again similar to that of 
Substituting Eqs. (3.27), (3.28) Into Eq. (3.22) and noting 
that (t) = -4>^ (-t) from Eqs, (2.30) and (3.23) we obtain 

(hi ) 

C2cosity-4Q,y(y+1) + 12QiY>3Q, + Q^] 

(2h^)^s1nTTY I I ^ 


’(h-j+x^)^ (h^-x^)^ 


■]= Fl(Xi) 


(3.29) 


where contains all the bounded terms. Now multiplying 
Eq. (3.29) by (h^+x-j)"'^ and letting x^-i'-h^ or multiplying by 
(h^-x^)''^ and letting x.|-»-h^, we obtain the following charac- 
teristic equation for the unknown constant y! 

2costiy + 4Q^(y-1)" - (Q-I+Q2) = 0 • (3.30) 
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One can see that this equation Is the same as those found 
In [5], [21]» [22], [23]. Eq. (3.30) is solved for Yi the 
power of singularity, which is assumed to be 0 <y<1 to get 
an integrable singularity. We can find this real value of 
Y numerically, for evsmple, using Newton's iteration method 
[24] within the desired accuracy. 


3.4.3 Solution IntaaKaZ Eouat-ion 


Define now the dimensionless variables <ii,T by 


= h-jO) , t = h^T , -h<|<x^,t<h^ . 


1 i 


(3.31 ) 


Then Eq. (3.22), with Eq. (3.23), takes the following form 
1 fl 1 4>l(h,T) 


p^(h^w) 


(3.32) 


Introducing the new unknown function 

4Pi <|)i(lif) 


in which p.j » p-j(x^) = constant, one can write 
fl , Oi(t) 

J Trh^k^^g(h^w,h^x) + iTh^k^^f(h^w,h^T)] dx = -v 


(3.33) 


(3.34) 
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Note that (1 -t^)"^ is the weight function of the Jacobi 
polynomials (t) . Thus the integral in Eq. (3.34) 

can be evaluated by using a Gauss-Jacobi integration for- 
mula given in [17]. Hence Eq. (3.34) becomes 

(j=l N-1) (3.35) 

in which are the roots of the Jacobi polynomials: 

(-Y.-Y) 


P„ (t^) » 0 . (i=l N) , 


(1-Y.l-Y) 

j 

N-l 


(3.36a.b) 


(o)j) = 0 , (j=l N-l) , 


(-Y.-Y) 


(Ti) 


and W^(i = l ... .,N) are the weights of Pj^ 

The singl e-val uedness condition* Eq. (2.21), which can be 
written as 

J^e,(T,)w, = 0 


(3.37) 


in this case, completes the system of (N) equations for (N) 
unknowns, 0^(T.j). We can replace this system by an (nxn) 
system by considering the symmetry of the strip about y- 
axis. Hence 

(j=l n) , (3.38) 


in which 2n = N and 





k„f(oij.T,) . l>,Ck„f(h,yj.h,T,) . k„f(h,a.j,-h,T,)] , 

, (3.39a,b) 


Is obtained. 

Examination of stress expressions Indicates that the 
stress components are bounded for this case except 
which has a singularl y at The definition of the 

stress intensity factor at the tip of the crack (x^«h-|, 
Xs^-hs) Is no longer the same as those In Eqs. (2.60), 
(2.73) since the power of the singularity Is not 1/2 any- 
more. Now define 

'‘a ■ 

which can be expressed as (Appendix E) 

ka “ -Q3(h^)^Pl0l(l) . a=h^ . (3.41) 

Qg being defined by (E.ll) In Appendix E. After solving 
the system given by Eq. (3.38) for 6 -|(t^), (1 = 1,. ..,n) we 
can compute 0^(1) using an appropriate extrapolation tech- 
nique. According to [18] 

m-1 (-y.-y) 

e,(l) = I c.P. (1) , (3.42) 

1 j=l 0 0 
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c. • i 

j RJ 1m j 

Val^Sfe 



(3.43a-c) 


Mm 


■ .■r-Y;-YP|f|pR.“y5 


m 


r"m+l 


'(tJ 


Hence 9^(1) can be computed after determining y and comput> 
Ing (1*1 ,,n). 

3.5 Stresses 


We consider the stress components which seem t.o be the 
most significant ones for a<h^ b=0. The stresses at 
Xg * -h^ can directly be obtained from Eqs, (2.90) as 


a2,,(-hi.y)/Pi * 


02yy(-^^yy)/P^ i 


- ^ 6l(Ti)h2i(p.Ti) , a<h^ , 
' 1 " 

[f j^e,(T,)h22(p.T,) . 


a<h 


1 ’ 


( 3. 44a-c ) 




T2,j,(-h,.y)/p, 


■ I ’ “‘''l ’ 

I 

■ "S' ’ 


where ( 1 = 1 , . . . ,n) are defined in Eq. (2.55a) for a<h^ and 
in Eq. (3.36a) for a = h^, From Eq. (3.8), by using Eqs. 

- 61 - 


(3.10) and (D.5), one can write 

4u, ra 


‘♦y, to 

> 2 yJ* 2 . 0 )*i^ k 2 l(» 2 .t)Sl(t)<lt . a<h, 

I ta 


(3.45) 


and following a procedure similar to the one followed In 
section 2.4 he further write 




^ « a<h^ • 


In" . 

I I^e^(Ti)Wik2^(n»Ti) , 


a»hi 


(3.46) 


where are the same as stated above and X 2 = h^n. 


From Eqs. (3.3) we obtain 
4y.> 1 (-a 
'l‘ -a 


a3yy(x3»0) = I k3^(x3,t)G^(t)dt , 


(3.47a,b) 




where 


and 


f" ^ -s(hi-t) 

k3i(x3,t) = J K3^(x3,t,s)e ds * 

-s(h,-t) 

h34(y,t) = J H34(y,t,s)e ds , 

KjiCxj.t.s) = I N,|^(s,t)[(a3|,+2C3|^ta,|,sx3)cosh(x3S) 

^ I 


(3.48a,b) 


+ (2a4|^+C4j^+C3|^sx3)s1nh(x3s)]e 
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I 


H3^(y,t,s) » i N^|^(s,t)[(a3|^-a^j^sh3)cosh(h3s) 

1 

-Shg 

- (c^|^-C 3 |^sh 3 )sinh(h 3 s)]e cos(ys) 


(3.49a.b) 


Defining X3=h35» and using Eqs. (2.49), (2.53) for a<h^ and 
Eqs. (3.23), (3.33) for a = h^ , then replacing Integrals by 
appropriate summations, we can rewrite Eqs. (3.47) as 


<>3yj,(X3.0)/p, 





a<h^ 


, a=h^ 
. a<h^ 


• a=h^ 


(3.50a,b) 


where are given again by Eq. (2.55a) for a<h^ , by Eq. 
(3.36a) for a = h^; r*hip, and 


l^3l(5fT^) ~ a[k3i (h35,ax . ) - k3i (h35»"a'f^)] p 


(3,51a,b) 


3.6 Numerical Results 

The problem Is solved for both plane strain and general- 
ized plane stress cases since the formulation is carried out 
with K, definition of which characterizes plane strain or 
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plane stress cases in itself (see section 2.1). Numerical 
results are obtained for material combinations II and III 
(see section 2.5 for definitions of material combinations). 
The cases where a<h^ and/or b<h^ are solved in the same way 
as the former problem. Solutions for these cases are ob- 
tained in order to establish limitations of the spring model 
approximation in the first problem. 

Figures 18-23 show comparison of the stress intensity 
factors for various geometries and material combinations 
for plane strain and plane stress cases. As one can see in 
Figures 18-23, the stress intensity factor does not depend 
on the material properties and plane strain or plane stress 
cases when h^'O. This is expected since h 3=0 corresponds to 
the case of colinear cracks in an infinite elastic solid. 
These figures show that the two solutions match very well 
for practical values of hj/h-j ratio (h^/h-j = 0-10%). Re- 
sults do not vary for plane strain or plane stress cases. 

But it is notable that agreement between the two solutions 
is better for material combination II. As h^-^® the stress 
intensity factors for a,b = 0,0.9h^ become the same for 
material combinations II and III, for plane strain and 
plane stress cases in the problem where the adhesive layers 
are approximated by distributed springs. Note that h^'^-® is 
equivalent to £ 2143-^0 which make the problem be the same 
as a center-notched strip having traction-free boundaries. 
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However, 1n the second problem the stress Intensity factors 
depend on the conditions of plane strain or plane stress and 
the material combination. In Figures 22,23 the stress in- 
tensity factors for a = b = 0.9h^ (for plane ' . -^ss case) 
first increase as h^/h^ increases (up to ".lO), then start 
decreasing with further increases in h^/h^ . This, we be- 
lieve, can be explained as follows: When h^ is too small 
the two cracks in adjacent strips are close to each other. 

As h^ increases the tips of these cracks will still be close 
to each other for small values of h^ whereas we will have 
the effect of a thin soft layer between main laminates. 

This will increase the stress intensity factors. But if h^ 
continues to increase, the cracks will loose the interacting 
effect of each other which will make the stress intensity 
factors decrease. As a result, for small values of h^. the 
effect of increasing the thickness of the softer layer is 
dominant*, for larger values of hg, effect of ceasing inter- 
action between two cracks dominates. 

Figures 24 and 25 show the variations of the normal 
stresses in x-direction on the boundaries of the adhesive 
layer. The /lormal stresses in x-direction in the adhesive 
layer and the adjacent laminate increase as the crack prop- 
agates. However, they can be reduced by increasing the 
thickness of the adhesive layer. In Figures 26 and 27 the 
normal stresses in y-direction on the line y=0 are shown. 
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These figures show that the cleavage stress ahead of a 
crack decreases as the th1cknesr> of the adhesive layer In- 
creases. This 1s very Important as far as crack arrest In 
composites Is concerned. One can see In Figure 26 that the 
cleavage stress at a location ^2 = -h-| can be reduced to 1/3 
of Its value in a homogeneous medium by bonding strips with 
adhesive layers of thickness hj = 5%h.| for a®0.9h^. Compar- 
ison of Figures 25 and 27 shows that there Is a possibility 
of delamlnatlon along the line ^ 3 ° when the crack 
approaches this boundary. The normal stress In x-directlon 
Is much higher than the normal stress In y-direction at the 
beginning of the adhesive layer. Therefore, depending on 
the strength of the bond between adhesive and laminates, a 
crack approaching the interface may propagate along the 
boundary (delamlnatlon) rather than going through Ihe ad- 
hesive layer. 

In order to solve the problem for a=h-| we should first 
determine the power of singularity, y» from Eq. (3.30). By 
using Newton's iteration method y is found for various 
cases as shown in Table 1(a). Figure 28 and Table 1(b) show 
the variation of y with X = p^/P 2 for plane stress and plane 
strain cases when v^=V 2 = 0.35. From these results one can 
observe that the power of singularity is higher for plane 
strain case. The difference is larger for small values of 
X and it vanishes as X becomes larger. One can also observe 
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that y“0.5 for both plane strain and plane stress cases 
when X=1 for which the composite medium is equivalent to 
a homogeneous medium. y is lower than 0.5 for X<1 and 
it is larger than 0.5 for X>1 . Figure 29 shows the varia- 
tion of the normalized stress intensity factor, K = kg/p^h|, 
with X on a logarithmic scale for h^ = .OShi* = V 3 = 0.35. 

For small values of X{<*-1) K is larger for plane stress 
case by an amount of more than 50%, then for larger values 
of X it is almost the same for both cases. The value of 
this normalized factor is much less for X>1 than its value 
for X<1 in both plane stress and plane strain cases, van- 
ishing as X->«, Variations of K for a = h^ with hg/h-j are 
given in Figures 30,31. Figure 30 shows the results for 
material combination II whereas Figure 31 does for combin- 
ation III. One can observe that (a) K is slightly larger 
for plane strain case, deviation being small {the same 
trend in Y)t (b) K is larger (-4 times) for material combin 
ation III, (c) K Increases with Increasing h^/h-j ratio. 

Figures 32-37 show some of che calculated results for 
stress components. In Figure 32, the cleavage stress 02 ^^ 
at Xg = -h^ , y=0 varies as the crack in main laminate prop- 
agates. The stress increases as the crack approaches the 
adhesive layer. The value of the cleavage stress when the 
crack touches the interface (a=h^) is unbounded if h 3=0 and 
it is bounded otherwise. Presumably this is the most 
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important point in this problem. The infinite cleavage 
stress (for h 2 = 0 ) at the beginning of the second strip can 
be made finite by considering even a very thin adhesive 
layer. This finite value decreases as the thickness of the 
adhesive layer increases. Figures 33 and 34 show the varia- 
tions of the stress components <^2yy ^2xy ^2 " ”^1 

line in y-direction. These stresses decrease as y increases 
and they are slightly larger for plane stress case, 
becomes smoother over 0<y<“ with increasing h 3 /h-| ratio 
(relaxation). In Figure 35, the variation of <^2yy 
Xg/h.! at y=0 is shown for material combination II. ^zyy 
again decreases with increasing h^/h^ ratio. Variation of 
the cleavage stress ^^yy adhesive layer as the crack 

approaches interface is shown in Figure 36. 
becomes unbounded when the crack touches the interface. 
Figure 37 shows the variation of «^ 3 yy with x^/h^ on y=0 line 
for a=h^ . 
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IV. CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 


The effect of the adhesive layers on crack propaga- 
tion in composites has been examined for both plane stress 
and plane strain cases. Composite medium consists of main 
load-carrying laminates and buffer strips bonded periodic- 
ally through thin adhesive layers. The adhesive layers 
have been approximated by distributed tension and shear 
springs. We solved the problem for the cases where the 
cracks were imbedded in laminates and/or buffer strips. 

The problem has been reduced to a system of singular in- 
tegral equations and this system has been replaced by a 
system of linear algebraic equations which has been solved 
numerically. The stress intensity factors and some stress 
components have been computed and presented in Figures 3- 
16. We saw that the stress intensity factors at the crack 
tips have increased by taking the adhesive layer into ac- 
count. They increase as the thickness of the adhesive 
layer increases (or as Young's and shear moduli of adhes- 
ive decrease). However, the presence of the adhesive layer 
relaxes the constraints at the interfaces so that the 
cleavage and shear stresses in the strip on the other side 
of the adhesive layer decrease. This is important from the 
point of view of crack arrest in composite materials. All 
the relevant special cases treated in literature can be 



recovered from this solution which has been observed to be 
in good agreement with them. 

Then we attempted to solve the case of broken lamin- 
ates. We observed that the spring model approximation is 
not mathematically suitable for this case. Hence we in- 
troduced the problem described in Chapter III in which the 
adhesive has been treated as an elastic continuum. We 
solved the case of imbedded cracks in order to establish 
limitations of the spring model approximation in the former 
problem. We concluded that the spring model is good enough 
within practical ranges. Approximation is not reliable for 
very large thicknesses of adhesive layer or for compara- 
tively weak adhesives. The case of broken laminates has 
been solved without any major difficulty. The character- 
istic equation to be solved for the power of singularity 
has come out to be the same as those obtained in [5], [21], 
[22], [23]. Stress intensity factors still increase as 
the thickness of the adhesive layer ‘ncreases. However, 
the cleavage stress at the closer edge of the second strip 
(first strip is broken) becomes finite in this case whereas 
it has been found infinite by ignoring the adhesive layer. 
Special cases can be recovered from this solution too. See 
detailed numerical results and conclusions in sections 2.5 
and 3.6. 
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The work presented here can be extended in several 
directions treating the adhesive as an elastic continuum: 

(i) The problem can be solved for non-isotropic 
(e.g., orthotropic) materials which will be more realistic. 

(ii) The strips bonded by adhesive can be of differ- 
ent material s. 

(iii) Finite number of strips can be studied. 

(iv) It may be useful to solve the problem for the 
cases where crack continues to propagate into the adhesive 
or along the interface (delamination) . 

However, one should keep in mind that all these prob- 
lems will require considerable amount of time and labor. 

We hope that the present work will contribute to future 
studies in this field. 


Combination 


Power of singularity, y 



Plane Strain 

Plane Stress 

II 

.8689 

.8658 

III 

.7183 

.7060 


(a) 


''rv-35 

Power of singularity, y 


Plane Strain 

Plane Stress 

.0 

.3203 

.2617 

.001 

.3207 

.2623 

.01 

.3240 

.2680 

.1 

.3538 

.3155 

1. 

.5000 

.5000 

10. 

.7551 

.7552 

100. 

.9160 

.9138 

1000. 

.9732 

.9724 


(b) 


Table 1. Values of the power of singularity, y, 
for various material combinations when 
a=h^ . 
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Figure 1. Geometry of the composite medium. 



2.2 ^ 


i 



-76- 


J 


Figure 3. The stress intensity factor kg vs. a/h? for the 
crack in main laminate (Combination I). 




Figure 4, The stress intensity factor Jck vs. b/h? for the 
crack in buffer strip (Combination I). 
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Figure 5. Variation of kg with respect to a/hi for b=0.3hi 
plane stress case (Combination I). 
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Figure 6. Variation of kb with respect to a/hl for b=0,3hi 
plane stress case (Combination I). 



ure 7. Variation of kg with respect to 
plane stress case (Combination 



Figure 8. Variation of with respect to b/h2 for 
plane stress case (Combination I), 


2.0 
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= 0.Sh 















Variation of k|j with respect to (log) h 3 /h 2 for 
b=0.8h2, plane stress case (Combination I). 





;e .3 

Og„(-h, ,y)/p, 


Distribution of the stress component o2yy at 
x 2 =-hi for h3=.05hi, plane stress case 
{Combination I). 
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} 



T2,,(-h,.y)/p, 


Figure 14. 


Distribution of the shear stress t2xv at X 2 =-hi 
for h3“.05hii plane stress case (Combination l) 



(-h,.0)/p, 



Figure 15. Variation of the cleavage stress cr2yy(~^l with 
h3/hi for plane stress case (Combination I). 













Figure 20. Comparison of two solutions for b=0.9hi, plane 
stress case (Combination III). 




kb/pjV/B" 
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Figure 21. Comparison of two solutions for b=D.9hi, plane 
strain case (Combination III). 













Figure 25. Distribution of the stress component ^2xx 
X 3 =-h 3 for h3=0.05h] , plane stress case 
(Combination II). 
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Figure 30. Variation of the stress intensity factor with 
h3/5»] for a=h) (Combination II). 








for plane stress case as the crack propagates 
(Combination II). 






Figure 34. Variation of the shear stress i2xy along the 
line X 2 =-hi for a=hi (Combination'^ 1 1 ) . 
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2.0 



Distribution of the stress component opvv at 
y=0 for a=hi (Combination II). 
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APPENDIX A 




Evaluation of some integrals [21], [25]: 

[^e"%os(ys)dy « , r>0 , 

»« 

J e“'^^sin(ys}dy « . »'>0 . 

|^ye*''-^cos(ys)dy = . r>0 , 

»w 

J ye“’’^s1n(ys)dy = -^pr|pjr . • 

/QO 

Jq rTrW)’ sin(ry)dr = (1-e , 

T yifen T sin(ry)dr = ye"®^ , 

lo VW ‘ ^ • 

#00 2 

Jo T?^JpF '“'‘('■J')'''' = Vs f*sy)e'*^ 

r cos(ry)dr = ^ (l-sy)e'®^ . 

0 ^ ' 


(A.la-d) 


(A.2a-e) 
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APPENDIX B 


i 


Expressions of the functions used in Eqs. (2.16) and (2.17); 

K^-3 sin(h4r) 

i4irHZs‘ + -^ (rHs^)] 

'0 

cos(h4r) 


4-2 cos(h4r) 

^ 3 (*> ■ VJ ^^ t * 

, Ac f“ sin(h4r) 

^4<®> = T TrTfiV • 


The forms of the functions a^j, b^j, Cy, d^j, 

( i=l ,2; j«l -4) appearing in Eqs. (2.17) are as follows; 


(B. 1 a>d) 


“ ^^ 32^35 ^31 ‘ 
^12 " ^^ 31^37 ■ ‘^ 5 ‘^ 32 ‘^ 38 ^^^ * 

“ ^^^ 32^39 " ^ 31 ^ 37 ^^^ * 
‘^14° ^^31^34 ■^‘^32^40^^° * 

^^21 ■ ^ V 30 ‘^ 35 '^ ^ 33 ^ 34 ^^° ’ 

^22 ■ ^^ 5 ^ 30^138 " ^^ 33 ^ 37 ^^^ ’ 
^^23 ^^ 33^37 " ‘^ 30 ^ 39 ^'^^ » 

*^24 "" “ ^^30^40'*' ‘^33^34^'^*^ * 
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“n ■ 'll ■ ■ -«,2 . 

‘*13 “ ■‘>32‘l42)/B . 

'*14 ■ <''3l''43*‘'32''44*/“ • 

'*21 ’ '21 > '*22 ■ -'22 ■ 

‘*23 ‘ ('>30‘'42*‘l33‘'4I>/“ • 

‘*24 ‘ - <‘I30‘'44*‘'33‘'43)''“ • 

*21 ' - {X,q27<:„ + q36'2i)/qg . 

'22 * ■ (»1<I27'12 * <l36'22>''‘)5 • 

423 = (X,q,- 11,4270,3. qjgCjjl/qj . 

'24 ‘ (^1<'2 - ^<l27'l4 ‘ ‘'36'24>''<'5 • 

*’21 ° '21 ■ *>22 ■ -422 • 

•>23 = - <‘'l*^‘'27‘'l3*‘’36‘'23>/‘>5 • 

I>24 = (^2 - ^l‘'27‘‘l4 ■ ‘>36‘’24)/45 . 

'11 = < M4 ' 21 -'•' i 1 i ': ii * X2 ''’ 2 ‘' 3 ' 21>''''2 • 

4,2“ (^2V22-4''l‘'l'l2*X2Sh2q3C22)q'’ . 

'13 “ '-'+W23-'''l‘>l<=I3*»2'''2V23>/‘l2 
'l4 ° *^2V24 -''’i‘'i'i4'^^2'''2''3'24>/‘'2 ’ 
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‘‘n * ®n » ‘*12 * "*i2 » 

' <^2 ♦ W23 ■ * XgShgqjdgjI/qg , 

^14* ^^2*'4^24"®^1‘*1**14'‘’^2®^2‘*3‘*24^^^2 » 
where 


° ■ Wai ■ ‘'32*133 ’ 

q, • (1-e‘^'''')/2 . qj . (lte‘^®''')/2 , 

qj • (l-e'^®'*2)/2 . - (l+e'^*''2)/2 , 

- 2 s(h^+h 2 ) -2s(hi+h«) 

^5 ■ El-e ]/2 . qg « [1+e M/2 . 

. -2sh9 "2shi -2shi "2shp 

q; » (e ^-e ‘)/2 . qg = (e Ue ^)/2 , 


^11 


-4shi 


q,- (l-e )/2 . q,;,. (1 h‘''*'’')/2 . 


( l - e -‘'%2 . q ,2 - ( He '’“"^)/2 . 


-4sh2, 


*^13 “ ^l+'^l)(X2-l)q7/2 


^14 * (^+K^)X2ShQ{X^qgq^-X3q^q3) . 

^15 * • 


^16 ” ■ ^4^^2^ * 

~ * *^18 ' ‘^2^5" ^7 ' 

q^g = -X4ShQq2q^ , qgQ = ' 


(B.2) 


(B.3) 
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-2shi 

- (l+K,)sh^qgC(X3-X4)e -{X3+X4)q^Q]/2 

*^24 * ®^2® * ^25 " (l"^^)shiqg « 

qgg • -(Ki+Xi)q5qg/2 . qg; • sh^e"^®*'Uqg/2 
qga - AXjX^s^^hJqg , 

‘’29 ’ ^X^sh^CXgq^qg+X^qgq^) . 

'‘30 “ ‘‘13 ^ “14 • 

-2sh2 

‘‘31 ■ ‘‘15 ■’■ ‘‘16® * ‘‘l7‘‘l8 '•■ ‘‘19‘‘20 * 

‘‘32 “ 0“^2^^®‘’2® ’ ‘‘22 **“ ‘‘23 " ^^2^3^45 = 

‘‘33 ^ ‘‘25® * ‘‘Zl " ‘‘27^‘‘28+‘‘29^ ■'■ ‘‘26 * 

% = (Xg-l)q4 - 2XgX3Shoq3 , 

‘‘35 ' ‘‘2 ■*■ ^^3®‘’o‘‘l ’ 

‘‘36 ' "‘’2“6 ^ t‘‘5"‘‘7^/^ • 

‘‘37 ■ ^^2“^ ^‘‘3 ■ ^^2^4^‘’o‘‘4 ’ 

‘‘38 “ ‘‘l ■*■ 2^43h^jq^ , 

‘‘39 = “l‘‘29 ^ ^X^sh^jqgqg + (Xl+4X3X4S^hJ)q^qg 


$ 


» 


^o‘‘24 » 
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^42 ■ *^1^*^5'*’^2*^29^ ' 

^43 ■ X2{[X2(l"^^3^4S*h*)“l ]q^ - » 

*144 ■ ^2^'^5'^^2‘^29^ ’ 


The expressions for the functions K^j(x^»t»s), 
(i»J“1i2) which appear In Eq. (2.33) are 

irK^^(x^,t,s) »^i^N^,^(s,t)[{a^l^+2c^(^)cosh(x^5) + c^^sx^s1nh(x^s)]e 

4 -shi 

TTK^2<^i»^*s) " I N2|j(s,t)[{b^l^+2d^,^)cosh{x^s) + d^|^sx^sinh(x^s)]e 

k*l 


I 


APPENDIX C 

Derivation of the stress intensity factors defined 
by Eqs. (2.60) and (2.73) as: 


^lyy(Xi,0) , a<h^ 

(2.60) 

^2yy(x2*0) , b<h2 . 

(2.73) 


From Eq. (2.31) we can write 


4Ui 1 GJt) 


(C.1) 


where ®jyyb(^j»0)» (j = li?-) are the bounded parts of the 
cleavage stresses, °jyy(^j*0)- Consider Eqs. (2.49) and 
(2.62): 

Define the sectionally holomorphic functions, 

•t rS. Gj(t) 

Xj(z) « 1 j J dt , (j=l,2) . (C.3) 

Following Muskhel ishvil i [12, Chapter 4] and using Eq. (C.2) 
we obtain 


1V2 

<J)4(-aj)e , 

X z = -J i + X. (z 


(C.4) 
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(C.5) 




/JaJ /SyiJ 




(J-1.2) 


where Xj^(z) and Xj^(xj) are bounded functions for 

Now with Eqs. (C.l), (C.3), (C.5), Eqs. (2.60) and (2.73) 

become 


4y, 

' T+icij" 


11m 

x^-^a 


4»T(a) 


4Mj 

-Th^ 


11m 

Vb 




a^hi 

b<h2 


(C.6a,b) 


Using Eq. (C.*^) we can rewrite Eqs. (C.6) as 

k lim /2{a-*,) ®1<*1> • 

® ' *'l x^-^a 

(C.7a,b) 

bb “ - ,y4,-'2lb-«2) S2(»2) . b<h2 

or with th : definitions (2.53) and (2.66) 


kg = - /a P] 0^ (1 ) , a<h^ 

kjj = - /F pg ^20 ) . *><^2 

are obtained. 


(C.Sa.b) 




APPENDIX D 


Expressions of the functions appearing In Eqs, (3.5) and 
(3.6): 




sln(h^r) 


r,j(s) 

F,3(s) 


TT 

4sf. 

tr 


r” K,+l 

cos(h,r) 


(r*+sM] 

dr , 


cos(h^r) 

r(r2+s2)* 


dr 


(D.la-d) 


^ 14 ( 3 ) 



s1n(h^r) 

(r2+S2)2 


dr , 


(1=1.2) 


.ynctions a^j, c^j, d^^, (i».j = l-4) used in Eqs. (3.6): 


®31 

= ‘'34/'>l 

, 832 » 

‘•35^“l ’ 

®33 

* >q29^°i 

. 83 ^ 

- Aq^j/D, . 

•>3J 

° ®3J ’ 

(j=l-4) 

9 

'31 

' ‘’32''“l 

’ ^32 

^33/“^l ' 

^33 

-Aq^g/D, 

* ^34 ‘ 


^3j 

° '3j ’ 

(j=l-4) 

9 

841 

° ‘'36/“2 

’ ®42 " 

^37^ ^2 ’ 

®43 

= -Aqj^/D, 

! ’ ®44 

= -Xq27/D2 
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''4j ’ -®4j • . 

*^41 “ ‘^33/^2 ’ *^42 * ^39^ °2 * 

*^43 * ^‘^30^*^2 ’ ^^44 “ ^**26'^®2 * 

^4j “ -‘=4j • • 

Cn = -(q2 + q4iC3^ + q42®4i)/q4o * 

^12 = (qi * ^41*^32 • ‘’42®42^/^40 ♦ 

*^13 “ ■ ■*■ ‘^41^33 ■*■ ^42®43^'^^40 » 

c^4 = - (Aqg + q4]C34 + q42^44^^% » 

'^Ij ' * ^‘>4l‘*3j ^ ‘^42^j)/'»40 * 

®n “ ^ 43 ‘^n V^31 ■*■ ‘^45®41 * 

^12 “ *^43^1 2 ■‘' ^^44*^32 ■*■ ^45®42 ’‘' l/t^”1)q2 
®13 ' ‘^43*^1 3 * %4^3Z %5^43 ~ ^/^^"'')q2 

®14 ■ ‘^43*^1 4 '*' *^44^34 ’’’ ‘^45^44 * 

hj " ^43^^^ ^ “ ‘’45®4j * (J=^-4) 

®2j "" • (-3=1-4) , 

'^2j “ ^ij ’ » 
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‘’2J ■ 'ij • 


(D.2) 


where 


■ 2(l-X)C(l-X)sh3e 


-2sh, 




•2sh 


+ (1+K^)A[(X-I)sh3qge ‘*-K3Xq5q7/2-K3q|,/2.q2/2] , 


-2sh, 


O2 • 2(A-l)C(l-X)sh3e '*-(l+AK3)q3q^]q2^ 


-2sh 


+ ( 1 +K^ ) A[ 0 -A )sh3qge -K3Aqgq7/2-iC3qf 2/2-q Jq/2 ] 


-2shi 

■e ')/2 , 

It 

CVJ 

or 

-2shi 

0+e )/2 , 

-2shq 

•e n/z , 

"4 = 

-2sh, 

0+e n/z , 

-4sh, 

•e ')/2 . 

^6 = 

-4sh, 

(1+e 0/2 , 

-4sh, 

e n/2 , 

^8 ' 

-4sh, 

0+e 0/2 . 


- 2 (h^+h 3 )s -2(h,+h-)s 

qg » n-e ^ M/2 . q^o = [1+e ^ ^ VZ , 


- 2 sh 3 -2h,s -2h,s -2h,s 

q^l = (e ^ -e ' )/2 . q^2 = C® +e ^ ]/2 , 


q^3 =“(l+K^)Aq2/2 , q^^ = (A-l)q,|^ , 

q^5 = (A-l)sh3q3 ^l’*’'^3)^q4/2 . 

qi5 = (1-A)sh3q^ - 0+K3)Aq3/2 , 

q^y = (1-A)q3 , q^g = (l+K^)Xq^/2 , 

q^g ® [1 - 0-K3)A/2]q3 + (1-X)sh3q4 , 
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(D.3) 


I 


^20 * (^-I)sh3q3 - [1 - (l-K3)V2]q^ • 

"2sh-i 

'^21 * * 
qg2 “ * 

qZ3 = -sh3qg - q^g/^ + ‘^3^12^^ * 

«»24 ’ ‘’14‘’21 ' ‘^9"’l3 ’ '’25 * '’l5'’21 ' '’l3'’22 » 

^26 “ '*16'’21 ■ '>13'’23 ‘ ^Z7 ' ‘’l7'’21 ^ ‘’l0'’l3 » 

*^28 “ '’17'’21 " '’9'’l8 ‘ '’29 " ‘’19^21 ' '’l8‘’22 ’ 

‘’30 = '»20‘’21 ■ '’18'’23 ’ ”31 " '’l4'’21 ^ '*10'’l8 ‘ 

932 ' n-^)^V 2 i * » 

qS3 = (1-X)>^q3q2i “ (l+Ki)^(qn+2^2‘’3^‘’l^^ » 

^34 = ^[(>^’’)sh3q3 + (1+«3)^V2^‘’21 

+ 0+K^)X[(q^2'^2Xq^q3)sh3 + 0+iC3)Xq^q4 + q2q33q2/^ » 

‘*35 ^ ^UX-l)sh3q4 - [l-(l-K3)X272]q3>q21 

- (l+K^)X[(q^2“2^‘’2‘’4^®’’3 ‘ ^’"'^3^^‘’2‘’3 ‘’2‘’3^‘>1 * 

‘’36 ' ^^^“’^‘’3‘’21 * n+'«l)^(qi2'^2Xq^q3)q2/2 . 

= X(X-l)q4q2^ - (l+Ki)X(qi2"2^q2‘’4^'’l^^ ’ 

‘’38 " ^n'’-^^^’^3‘’4 ■ 0+'^3^^‘’3/^^‘’21 

+ (l+K^)X[Iq^^-2Xq^q4)sh3 - (l+K3)Xq^q3-q2q4]q2/^ ’ 
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I 


^39 • A{(l-X)sh3q3 ^ n-0-iC3)A/2]q^}q2^ 

+ 0+K^)X[-(q^^+2Xq2q3)sh3 - (l-K3)Xq2q4+q2q4]q^/2 , 

-2sh, 

^40 * \l”X)sh^0 * 

%] ’ (l+«3)qn/2 , q42 - n+K3)q^2^2 . 

q 43 - -sh^tanh(sh^) + (1 +k,)/2(X-1) , 

^^44 ’ -^^+*^3)q4/2(^-l)q2 . 

*^45 * 0+>«3)q3/2(X-l )qg . (D.4) 

The functions K^j(x^,t,s), (i»j®l,2) appearing in Eq. (3.10): 

^ *sli 

‘>TK^T(x^,t,s) = ^i^N^j^{s,t)C(a^l^+2c^l^)cosh(x^s) + c^l^sx^sinh(x^s)]e ^ , 

irK^2(>‘l '^*s) = I N^I^(s.t)C(b^j^+2d^j^)cosh(x^s) + d^l^sx^.sinh(xjs)]e ^ , 

K“ I 

(1=1,2) . (D.5) 
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APPENDIX E 


Derivation of the stress intensity factor for a*h-| , 
b=0 in section 3.4.3: 

Define 

Y 

kg = lim (xg+hj) O3yy(x3»0) , a»h^ . (3.40) 

’‘ 3'*’"^3 

From Eq. (3.47a) we can write 


°3yy^’‘3*°^ " k 1 | * '^30b^^3^ ’ (E.l) 

• -hi 


for a = h^ where cr30b^^3^ bounded and 

r -s(hi-t) 


(E.2) 


in which 


^31s^^3*^*^^ “ 1 im K 3 i(x 3 tt»s) 


S-K» 


From Eq. (3.49) by using Eqs. (2.27) and (D.2) 


(E.3) 


l- 2 (h^-t)s 3-2(h3+X3)s -s(h 3 +X 3 ) 


'^315^^3*^’®^ ^ ■ 2 X+K^ * 1+Xk 


]e 


(E.4) 


can be written with which (E.2) gives 


where 
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I 


Q5 ■ 2[{A-^»c^) - (I+Akj) ] 

•1 “1 
Q5 ■ C(X+K^) - 3O+X1C3) ] 


(E.6) 


Substituting (E.S) Into (E.1) we obtain 


3yy^^3’®^ ■ " 2u^ “ I CQgChj+Xg) Qg]Ct-(h<j+h3+X3)] G(t)dt 

3 


* '^30b^’‘3^ • 


(E.7) 


The Integral In the last expression can be evaluated by 
following the procedure followed In section 3.4.2. Hence 
one can get 


^3yy^^3»0) 


-2Q, 


“I*] (h] ) 

(2h,)Y 


(ha+Xg)'^ 




where 


(E.8) 


Q4 = Pi(Q5Y-Qg)/s1n7rY 


(E.9) 


and c^ 3 i(,(x 3 ) Is again bounded. 

Now substituting (E.S) Into Eq. (3.40) 
<l>^(h^) 


kg = - 2Q4 ^ 


(2h^)^ 


= -2Q. Htn }/7 (h,-x,)\(xJ 
^ I I I I 

= - Q 3 ( hi )^ PTe ,( l ) 


(E.IO) 
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Is obtained. Here 



04(2) (K] + U/m^ 


(E.ll) 
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% 


f 


